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ABSTRACT 


This thesis investigates stable, asymptotically stable 
and unstable dynamical systems from a topological point of 
View With direct application and interpretation to systems 
of differential equations. Knowledge of topology is not 
a prerequisite. Specifically, such concepts as Poisson 
and Liapunov stability as well as parallelizable and 


dispersive systems are investigated. 
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INTRODUCTION 


Phonic utowecurrenetyeavallable, stability theory of 
me Meonledis syauciiauregulres an extensive background in higher 
mathematics. The purpose of this thesis then is to make 
this theory accessible to the reader who does not have this 
Eer Eronda lulls asoumed that the reader has taken courses 
in elementary differential equations, intermediate mathe- 
farnedt analysis, and @amear algebra. All the topological 
properties needed are described in Section l. 

sections 1 - 3 deal with developing basic concepts of 
dynamical systems. Sections 4 - 6 use these concepts in 
studying various kinds of stable or unstable systems. Topics 
of recursiveness and dispersiveness are also discussed. 

Among the stable systems, Poisson and Liapunov stability 
will be discussed as well as non-wandering systems. Poisson 
instability, Lagrange instability, and complete instability 
will also be investigated. Topics in dispersiveness include 
parallelizable and dispersive systems. The relationship 
between several of these systems will be shown. One such 
relationship to be shown is that a parallelizable system is 
dispersive, completely unstable, Lagrange unstable and 
Poisson unstable. 

several important and revealing theorems will be pre- 
sented. Among these are the theorems involving the Liapunov 
function which allows analysis of a system as to stability 


or instability without actually solving the system. 


L 





The book Stability Theory of Dynamical Systems by 
nata and SZego Wilibpe used extensively throughout this 


thesis and 1ts bibliography contains an excellent list of 


references. 





IA EEES TCAL- PROPERTIES 


DM ero o als section is to expose the reader to 
all of the topological properties needed to understand the 
material presented in this thesis. It is important that 
the reader develop a firm understanding and familiarity 
with these concepts to facilitate in understanding of the 
material to follow in later sections. 

nn zreeinamıch chezesncept of distance. 

Po Ss AMO Emp set, a metric Or distance function, 
on X is a real-valued function d of ordered pairs of elements 
Sica sablisty theeesconditions. First, d(x,y) >0, 
Suomi x.y) = Of andeonly if x= y. Second, d(x,y) = d(y,x). 
foes Called the symmetric property. Finally, 

Mere aa) + d(Z.y). This property is called the 
ete inequality, “A jmeuric space X is then a set X with 
a-metric defined on it. 

If Xo is any point of a metric space and r is a positive 
zeal number, the open sphere S(x,,r) with center x, and 
radius r is the subset of X defined by S(X_ sr) = {x:d(x,x,)<r}. 
NS US o As seta to be open if and only if it is 
the union of open spheres. Since x, is clearly an element 


0 
of S(xy>T); we say that S(X,>T) is a (spherical) 


Meteoro son XQ: 
If A is a subset of a metric space X, a point x in X is 


Cr Commerce Olillauylon point of A if each open sphere 





centered on x contains at least one point of A different 
from x. The essential idea here is that the points of A 
different from x get arbitrarily close to x. Then a subset 
Fof a metric space X is called a closed set if it contains 
are swaree mussen Ponts Another characterization of 
lose Set is thatea seu E G X is closed if and only if 
Meow complemeny ai swan Open set. Then by definition, 
the whole space X and the empty set Ø are open sets. Then 
Baches emptey set isthe complement of the whole space, 
we have by our previous observation that the whole space X 
and the empty set Y are also closed sets. 

Now, in any metric space, the union of any collection 
of open sets is itself an open set and the finite inter- 
error of open setszıszan, open set. Taking set complements 
it follows that any intersection of closed sets is a closed 
a FE mestcheseinıiverungensei Closed Sets is a closed set. 

TP Xo is a point in the metric space X and’r is a non- 
negative real number, then the closed sphere S[xy>r] with 


eenter x, and radius r ís the subset of X defined by 


0 
S[xy,r] = tx:d(x,x,) < r}. 

Let X be an arbitrary metric space and let A be a subset 
Ci a Called an interior point of A if it 
is the center of some open sphere contained in A; and the 


interior of A, denoted by int A, is the set of all its 


X 
interior points. Symbolically, we have 
intyA = ads.) e A for some r}. Since the 


interior of a set is a union of open spheres, it follows 
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that intyA is an open set. Hence we have the characterization 
Mata set Ads Open if and only if A = intyA. Wezor ven 
write int A an place of intyA whenever no confusion will 
ee sult . 
Et X ismaremetric space and A is a subset of X, the 
Gesture Ole am A, denoted Dy clyA, Mane Union of A and 
ae set oral its accumulation points. Intuitively, clyA 
IS Aitselt together with all other points in X arbitrarily 
akose vo A. Since it as clear that cl A Contains all its 


ecemmulation points, cl A is a closed set. Thus we nave 
an additional characterization of a closed set: AGC X is 
emoscd if and only UA = cl yA. We often write cl A in 


place of cl A. 
Again let X be a metric space and let A be a subset 
of X. We will use the notation X-A to denote the complement 
AO IA DONT Ins called a boundary point oí A 
if each open sphere centered on this point intersects both 
A and its complement X-A. Now the boundary of A, denoted 
bays. 1S the set of ali@the boundary points of A. It is 
Pret ly shown that bdy A = cl,A A cl (X-A). Moreover, the 
boundary of A is a closed set, and a set AC X is closed 
Everyametriesspace has, in addition to the properties 
already stated, two important "separation" properties. 


First the Hausdorff property: aif x and y are two distinct 


Ponts In) Xk, vem there exist disjoint open sets un and oe 





MeN thay xX £ ur and y € U, Second, the normal space 
Property: 1% A and B are disjoint closed sets of X, then 


there exist disjoint open sets U, and U, such that A CU 


A B A 


and BE Ur: 

We now use these basic definitions to describe more 
advanced ideas in topology. We are interested in character- 
izing "compactness", "completeness" and "connectedness". 

ec be ally metric space. A class {G,} of open subsets 
Nos adobe atimepen Cover of X if each point in X 
belongs to at least one G,, that is U G, = X. A subclass 
of an open cover which ís itself an ee cover is called a 
pupeover. Mhen a compact Space is defined to be a metric 
space in which every open cover has a finite subcover. This 


becomes a familiar definition when one recalls the Heine- 


Bewer beorem umich states that if [a,b] is a closed interval 


[a,b], then there exists a finite number of the G,'s whose 
Bon contains [a,b]. This is just the statement that 
closed intervals Sn the real line are compact sets. 

A Mur Deus oO prepersy Of a compact space is that any 
closed subspace of a compact space is itself compact. Another 
useful property is that any infinite sequence in a compact 
set contains a subsequence which converges to some point 
ieee Comiacemsey, this also becomes a familiar property 
when one recalls the Bolzano-Weierstrass Theorem: Every 


era <eumalesron. zzrnis theorem is just an application of 





this property applied to the real line. These last two 
properties are very important and will be used extensively 
eozoushout this thesis. 

To define the notion of a complete metric space we need 
the concept of convergence. If iS a sequence in X, 
we say that C Icon VE Sent vO the point x In X if given 
ay pe =) 0, there exists a positive integer N such that 
d(x, 5x) < £ whenever n > N. Then in particular we have 
d(x, 5x) < € whenever n,m > N. A sequence with this latter 
Property is ealleda a Cauchy sequence. Now, not every Cauchy 
sequence is a convergent sequence. For example, consider 
Mmeewio! Towine Sequence an the metric space of rational 
mime cum see os l, selene s-241, 3.1415, ...}. Clearly, this 
tomas cauchy sequence but IT is not convergent because it 
SONVerges CO 1; DUb T is not a rational number. Motivated 
by this example we say that a metric space is complete if 
every Cauchy sequence in X is convergent. Completeness is 
related to compactness in that a set is compact if and only 


A e A ee 


eat Tals complete :and wVowally bounded. By totally bounded 

we mean that the space X is the union of a finite number of 

Meee ovmercs Of radius less than r for each positive r. 
Finally, we define the concept of connectedness. A 


Eorneceredespeeris a Metric Space X which cannot be repre- 


sented as the union of two disjoint non-empty open sets. 


mos AU B where A and B are disjoint non-empty open sets, 
then X-A = B and X-B = A so A and B are also closed sets. 


Isa Spertci s connected if and only if the only sets 


O E 








Sms, ser P2 This property is very important and it can be 
ONO thar the intermediate value theorem in calculus is 
equivalent to this property for intervals on the real line. 

WEN NOW tateention to the concept of a continuous 
function. This concept has several different characteriza- 
rons and wenvepine withthe basic definition that looks like 
the definition used in elementary calculus. 

fom ce oetuneuemom trom a metric space X with metric d 
to a metric space Y with metric d'. Then f is said to be 


S ae a e or each e > 0, there exists a 


0 
ô > 0 such that for x e X, we have d' (f(x) ,f£(x_)) ne 
whenever d(x, xq) < 6. The function f is said to be continuous 
As CONTINUOUS at cach point of X. The definition we 
Tun calculus was Hor real-valued functions so the metrics 
just involved the usual absolute value for the real line. 
one char ererizations of a continuous function 
foe Demncededsin thismemesis and are stated here. A function 
REN NINAS anggonv it the inverse image of each open 
(closed) set is an open (closed) set. Also a function f is 
continuous if and only if for each sequence {x,} which 
converges to x, the sequence {f(x _)} converges to f(x). 


We now state some important properties of continuous 


functions 


A Continuous real-valued function f 





following sense: if a = inf{f(x):x e A} 


and b sup{f(x):x e A}, then there exists 


X} >X, € A such that f(x) = a, and 
£(x,) 


Dr 


The continuous image of a compact (connected) 


set is compact (connected). 


Next, we define the concept of a homeomorphism. A 
homeomorphism is a one-to-one continuous function from one 
metric space onto another metric space such that the inverse 
of this funetion 1s also continuous. (A one-to-one function 
es runerion which assıens different images to different 
ements in its domain. A function f from a metric space 
meena metric space Y ds onto if the image of X under f is 
the entire space Y. If a function is one-to-one and onto, 
then its inverse is well defined.) 

We have already stated the definition of a compact space. 
A closely vteleabted 1deau@secthat of a locally compact space. 

AS iS name suggests, a metric space X is locally compact 
if and only if each point of X has at least one compact 
neighborhood. Hence, a compact space is automatically 
Totally compact. 

One of the most important properties of a locally compact 
space is that through the addition of one point, usually 
PLeacomuetmuOomocmmneos1 deal point. the space can be extended 


to a compact space. This extension is known as the one-point 
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compact i eati on = ine Cne—-point compactification of a 

metric space X is the space X' = X U {w} where the open sets 
in X' are all the open sets in X together with all subsets 

U of X' such that X'-U is a compact subset of X. In other 
words, a set U is open in X' if and only if UN X is open 

and whenever w e U, X-U is compact. As this is an important 
fact, we give a formal proof. We must show that X' is compact 
and that unions and finite intersections of open sets are 


open. 


THEOREM: The one -porne compactification X' of a metric 


space A ls compact. 


ARSS. insbe intersections and arbitrary unions of open 
sets in X' intersect X in open sets. If wis a member of 


miemunyersection of two open subsets of X', then the comple- 


ment of the intersection is the union of two closed compact 


Me CS Ofk and is therefore closed and cómpact.  Hence 

the intersection of two open subsets is open in X'. If w 
belongs to the union of the members of a family of open 
subsets of X', then w belongs to some member U of the family. 
Moreover the complement of the union is a closed subset of 
the compact set X-U and is therefore closed and compact. 
Hence the union of open sets is open in X'. Consequently, 
arbitrary unions and finite intersections of open sets in 

X' are open. Next, let U be an open covering of X'. Then 


wis a member of some U in ù and X-U is compact. Hence 
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there is a finite subcover of U which covers X-U and thus 

this finite subcover together with U is a finite subcover 

cme’. Thus X! is compact. I 
Jur fina EOP an tnis section of topological prereq- 

tes SINE concept ol product spaces. A finite product 

Smece 1S the cervesian proquet ol two metric spaces.  1n 

ties thesis, one of the metric spaces will usually be the 

real line R; consequently if X is a metric space, we have 

the product space to be XxR. Hence any point in the product 

space will be written (x,t) where x e X and t e R. Now the 

ọpen sets in the product space are just the cartesian products 

of open sets in the metric spaces. That is if U is an open 

set in X anar is saneepemeses in R, Chen UxV is an open set 

in AxR, and every open set in XxR is of that form. If dis 

che metric for X, we define the metric d' on XxR in the 

following manner. If x,y e XxR, then there are elements 

a,b e X and r,s e R such that x = (a,r) and y = (b,s). Then 

Ber alas) 3 nz. The spaces X and R are called 

the coordinate spaces of the product and the functions Py 

and Pps which carry the point (x,t) of XxR into x and t 

ie Peoria yma re  Caibedecne projections into the coordinate 

spaces. It is easily shown that the projection mappings 

ASCO IAN SA Ne Cons and that their inverses are also 

continuous functions. Moreover f:Y + X is continuous if 


f and P of is 


and only if each of the compositions Py y 


continuous, Lor any metric space Y. 
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A function f = f(x,t) defined on the product space is 
San OUS nnd only if for cach t the function f( 
whose value at x is f(x,t), is continuous. Similarly, 
Meet) 1s Convimuous in t if and only if for each x € X, 
me function £(x. ) such that f(x, )(t) = f(x,t), is con- 


Anous. Finally if f is continuous on the product space 


then f is continuous in both x and t by our previous remarks. 


However, the converse is not, in general, true. 


15 





II. DYNAMICAL SYSTEMS 


We seek to motivate the concept of dynamical systems 
by some examples of differential equations. But first we 
need the existence and uniqueness theorem associated with 
such systems. 

We consider systems of first-order differential equations 


written in the form 


jr: 
dt = f(E; Xetea Xn) 
AX, 
a L (E5X 3. ..»Xx,). 


Such systems frequently arise out of physical problems where 
the independent variable t represents time. Then a solution 
of the above system describes the "trajectory" of a particle 
moving in n-space. To see this more clearly, we look at 

the following two examples. 


Consider the following system defined in R?: 


ax 
rn : 
ÓN -sint 
ax 
2R- 
at cost 


16 





a 


ee = 


q 
at 


Moen the general solution to this system is given by 


Xx = COS tet Cy 

X, = sint + Co 
= + 

X3 t ee 


where the c's are all constants of integration. 
Suppose we specify that this solution passes through the 


=a Ty X, = 0, and X3 = 0 at time t = 0. Then the 


Point X] 


solution is given by 


x] = cost 
X, = sint 
Xa = t. 

3 l 


Such a solution describes the "trajectory" sketched below. 
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It is clear that at any time, we may find the exact 
Per er son of the particle. mihe same is true of the following 


example. 
> 


(nz ıder thezsystemsderined in KR given by 


o 
ab 
dx 

2 
a O 
one 
at 


Mocre we Specie Chat its solution passes through the point 
Xy = 0; X, = 0, and x2 = O at time t = 0, Then this system 


has the solution given by 


18 





DD 


X, = t 
—— 
3 3 


Its "trajectory" is sketched below. 


K3 


E) 


Mancan na e a aae orob lem Tor such a system consists 


Ame electing an initial time t = t, and a point X, in Re 


O 0 


and Chen seeking a solution whose graph passes through Xo 


at time t = Go: 


By setting X(t) = (x, (t),...,x,(t)), 


A). O) 


this system may be rewritten as 
dX 


= f(t,X), where F (x, '(t),...,x,'(t)). 


ax 


We next want to examine conditions under which unique 


solutions ekist for an initial value problem of the system. 
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Suppose that Mit) 1s bounded and continuous in a region 
Dot ee Busen znoge, assume that í satisfies a Lipschitz 
condition in D; that is there exists a non-negative real 
number b such that whenever (t ,X,) and (t ,X,) are any two 
Besnts in R, d(f(t,X,), f(t,X,)) < ba(X, ‚X,) where d is 


the ordinary Euclidean metric. Then we have the following 


theoren. 


Boa ueber continuous, bounded, and 
Sawisivea D ehiz condition in a region D 


a 


of RUS, and let (t,,X,) be any point in D. 


Inen he mra value problenm 


=> = f(t,X) 


X(t) = Xo 


Do OUEN Oa Oon on an interval I of R 


Gr Se H 


ÁE o o ë: ë mmm AA i SS 


We will be interested mainly in systems defined on all 
or R. -Tnen by ereuriıng that’! is continuous and bounded on 
all of ne auossabisties a global Lipschitz condition, we 
will have existence and uniqueness of solutions to the above 
system. 

One elegant and simple proof of the above theorem requires 
a knowledge of contraction mappings and some theory of 


fixed points and will not be presented here [Kreider, Kuller, 


Ostberg, p. 385]. 
20 





oia sa coma mool may De found in Theory of Ordinary 
Differential Equations by Coddington and Levinson. 

Of critical importance is the question: what happens 
to the solution if small changes are made in the function f 
Dein the initial conditions Eo and Xp? In most physical 
problems whose solutions are obtained from such a system, 
experimental error and estimations of physical constants 
Mayeattect the initial conditions, and the presence of 
variable parameters, such as temperature or density, may 
affect f. We would hope that such small variations would 
Moumealisce a drastic Variation in the solution. It can be 
shown that these unique solutions guaranteed by the theorem 
dependa continuously on f and on t, and X, for a certain 


0 0 


enaunt of time; that is small changes in f or in Eo and Xo 


produce small changes in the solution. More precisely, 


Lome — 1 Cue emecmcomcinuous, bounded, and 
satlisty a ea condition in a domain D 


of the (n + 1)-dimensional (t,X) space, and 
interval I: a <t <b. There exists a 6 > 0 
such that for any (t, Xp) where a < t} < b, 
and d(X,,X(t,)) < 6 there exists a unique 
solmwson Yon I with Y(t,,t]>%,) = X, where 


a < t, < b. Moreover, Y is continuous Ina. 


(Proofs of the continuity of solutions may be found in 
Kreider, Kuller, Ostberg, p. 390 and Coddington, Levinson, 


Dee) 
an 





What happens to these solutions over an extended period 
Oc ime involves the theory of stability which is the primary 
merus Of this thesis. 

We will be interested mainly in first-order autonomous 
systems, that is, "time independent" systems. Such a system 


may be written in the form 


X= fee) where X denotes = 


X(t) = Xo: 
The following examples provide some insight into such 


autonomous systems. 


EXAMPLE 2.1 Consider a simple pendulum of mass m swinging 
w out friction on a weightless rod of length 1. -By 


Newton's second law, the equation of motion of the pendulum 


ls given by 
2 
-mgsin 9 = ml =- 
av 
or 
( 
6 = -g/1 sin 0 
2 es m 
By setting g/l = k” and 0 = w we 





obtain the following system: 
Se quilibrivm 
position 
ANS 
22 





w = 21° Sin 9% 


It can be shown [Kreider , Kuller, Ostberg, p. 397] that 
Be general solution to this system is 


we = 2K° COS O + Cc 


Where c > ap is a constant depending on the initial 
conditions: 


Graphing several of these solutions in the $w-plane we 


have the following diagram. 
6) 





PTCURE 2.1 


To see how this diagram was obtained, we see that w is 


the angular velocity. Suppose c = -2k°, Then when 8 = 0, 


25) 





w = O and hence 6 and o are both zero. Thus the system is 
Bor chaneing with time. The point i) is the complete 
meeTeeuory Tor this case. Physically this shows that if 
the mass is started at the equilibrium position with no 
Seomlar velocity, the mass will remain at the equilibrium 
Pes ction indefinitely. Suppose c = oyó and the system is 
started at 0 = 7. Then again w = 0 and hence 8 and w are 
both unchanging with time and the point > describes the 
complete trajectory. Physically this means that if the 
mass is balanced vertically above the equilibrium point 
MO angular velocity, it will also remain in that position 
indefinitely. Clearly the same analysis can be made for 
Rn NCh is exactly what we would expect. If this 
EN stem 1s started at any other angle €, it will have an 
taeerat velocity which will cause the mass to approach the 
angle 8 =7T or 8 = -—-T and it will get arbitrarily close 
ent will newer actually reach that point. Neither will this 
mass swing back away from this point but rather it will 
continue to approach the angle 0 = mt or 0 = -7 as time goes 
Panic This is the physical meaning of the trajectories 
between the points (¿”) and Cor 

Now suppose -2K° Sue DK, From the previous remarks 
it is clear that the mass will have a velocity at @ = 0 but 
Tpomveleciuymiscenol ereat enough to permit it to approach 
the angle 8 = T or 8 = -1T. We would expect in this case 
that the mass would swing back and forth across the equili- 


brium point always reaching the same angle on either side 
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of the equilibrium point. This is the physical interpre- 
Pan ron ot tne elliptical trajectories shown in the diagram. 

minallvy, suppose € > 2K. Then o is always positive 
which means that w is always positive or always negative. 
If w is always positive, then 6 is always positive and 
hence 8 is always increasing. A similar statement follows 
if w is always negative. Physically this is the odie uae 
the mass has been given enough initial velocity to swing 
completely over the top position. In a frictionless envir- 
emment we would expect this behavior to continue indefinitely. 
Bes che physical interpretation of the top and bottom 
trajectories in the diagram. 

To see the direction Of a particle along any trajectory, 
pick an angle 6. For 0 < 0 < 71, we see that o Oo commu as 
decreasing. For -īm < 0 < 0, we see that w >= 0 so Wis 
increasing. This completes our analysis of the diagram. 

It is noted that in this example the solution is not 
given as a function of time. To obtain such would mean the 
difficult task of solving a non-linear differential equation. 
Inis is not any more informative since the movement of a 
point along any trajectory with increasing time can easily 
be determined from the system describing the motion as we 
have just observed. 

A more elementary approach to the problem is by observing 
that for small angles 0, sin 8 is approximated by 8. Thus 


ou Eey Stem would be given Dy 
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It can be shown that the solution to this system is 


cD 
H 


Asin kt + Bcos kt 


Akcos kt - Bksin kt 


e 
ii 


where A and B are constants depending on the initial conditions 
[Spiegel, p. 210]. It is seen that this solution describes 
simple harmonic motion and closely approximates our trajec- 


tories for small c where _2x° oz 2K°, 


EXAMPLE 2.2 We consider the following vibrating spring 


problem. 


Equilibrium 
Ost el Owe le a 





The equation of motion is given by: 


SL + by + ky = 0 
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where W is the weight of the mass, g is the force of gravity, 
sc he damping constant and k is the proportionality 
constant of the spring. Suppose W = 6 lbs, g = 32 lbs, 
Nana ID 0 is left as an undetermined constant for 

now. 

Then it can be shown [Spiegel, p. 195] that the general 


lution for y is 
-Spt 8 8 | == 3 
y=e (Acos (= 19 nC eat Sit (3 O 


where the constants A and B depend on the initial conditions. 
iaewesspecity the inátial conditions to be y(0) = 1/3 


and y(0) = 0, we have the solution, 


8 
-—bt 
y = OL 3 (3 cos (S49 - pS t) t 


b sin(S 9-b? t)) 
os 


To make this example more explicit we will consider two 
Cases for b. First an undamped system with b = O and then 


a damped system with b = 1.5. With b = 0, we have 


y = cos 8t 


wir 


With b = 1.5, we have 
yes g € (3cos 4 f3 t + {3 Sin UN? E) 


By use of a trigonometric identity, this last equation may 


an 





be written 





y = 213 ort sin (4S3 t+ 2) 


Now that we have the solutions, we change this second- 
order problem into a system of first-order equations. Let 


X = > and y = x. Then forb= 0 


x(0) = (1,3) 


Mor notation, let the initial point TA = X» Then the 


solution to this system is 


-Ssin 8t 
X(X_ st) = ( y 
3008 8t | 
For b = 1.5, we have 
-bex _ key 
Mm = L(X) = W W 
x 


x(0) = (7/3) 


and the solution to this system is 
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os (Uv) p ze” "sincu St) 


Ne 


X(X >t) = 


Be "sinch St E 3) 


It is clear that as t + œ, X(Xp ot) > o which we would 
intuitively expect from a damped spring. 


We sketch these solutions on the x,y-plane. 





2 a 


For the undamped problem we see that its trajectory 
descrives simple harmonic motion. It continues to oscillate 
indefinitely with the same maximum displacement from the 
SONES cion. On the other hand, the displacement 
and its rate of change both decrease with increasing time 
for the damped problem. 


since we have two different systems, it is necessary to 


eg 





Moc their solutions on two different planes. This helps 


Bosawoid any misunderstanding in our future work. 


EXAMPLE 2.3 Consider the second-order differential equation 
fey - Of, where y is considered as a function of time. 
settings X = y, this equation may be converted to the 


system 


Now considering X = (2) Marne iia condition that at 
NU Ne SOLUTION passes through the point fk we have 


the following initial value problem 


x = f(x) 


ly. 


X(0) = Go 


where f is the linear function defined by 


f,(t,x,y) 
a) = e (t,x,y) = 


EDF A simple check shows that 


this system has the solution 


E OS 
X(t) = cam )° 


If we let X, denote the point > then X, denotes the 
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vial conaltion. Further since the solution depends on 


this initial condition, the solution X(t) is more properly 


cos a 


stated X(Xy;t) = Coin + 


There are three observations we wish to make about this 


solution. First, we see that X(X_ 50) = Xə} Next, 


X(X(X_,t5),t,) = X(X_ st = t,). 


1 


cos Us 
This is true since X(X(Xp,€62),t,) = X(( 


in ur,“ 


The term on the right is the solution of the system which 
COSME 


0 = oe E TU=*So lion is 


Mase the initial condition Y 


cos to cos t - sin to sin t 


X(Y¿>t) = ( ). Hence 


cos to sin oet sin to COSAC 


cos Us COs Ez - sin to Sib Et 


Yo) = Coo, t, sin t, + sin t, cos er Using a 


trigonometric identity we have 


cos (t, + =. 
) e 


2 Thus a 
sin (t + to 


K(X, ,t,) = ( 


AO tdt IS X(Y,,t,) = X(X).t, + to) 


is continuous because both of its components are continuous. 


Finally X(X_ st) 


sie that the graph of the solution is a circle 
in the x,y-plane, and if we consider the graph to be the 


EuaceuOryeorechie point X that point moves counter-clockwise 


0? 
with increasing time. This can be seen by choosing some 


positive value for x. Then y is positive so y is increasing. 


Sul 





RES 


C 


Eee aler the initial condition to X, > b 


0 ), the 


a cos t - b sin t 


lution is X(Xp t) = ( Wand ibs trajectory 


ASE F D COS L 


sta circle passing through the point Xo’ 

Now referring to Example 2.3 we see that the solution 
X(X,>t) is a map from the product space RÓXR into the space 
E. We note also that Ro Sima me ias paces. Lhe three 
observations made about the solution are precisely the 
Meeoperties defining a dynamical system. The fact that the 
system in Example 2.3 does define a dynamical system will 
Bemimade clear by the following definition where in this case, 
the map m is defined by (x,t) = X(x,t) where x denotes a 


Pont in aoe in all further work, X denotes an arbitrary 


metric space. 


DEFINITION 2.4 A dynamical system on X is the triplet 
(X,R,T), where T is a map from the product space XXR into 
eR pace X Satisfying the following axioms: 


oe very x In X, T(X,0) = x, 
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A E SE ia and t-ta in- R, 


aa? 


m(m(x,t,),to) = m(x,t + to), and 


1 

3. The map t is continuous. 
Hee space X and the map TT are called the phase space and 
mee whase Map, respectively. 

Demeter szpeint on, whenever no confusion will result, 

the symbol t will be suppressed. Thus the image under t 
of the point (x,t) will be written more simply as xt. The 
first two axioms then become 

IO every x in X, xO = x, 


2. For every x in X and t nR, 


122 


xt, (t,) = x(t, + ts). 


ALIS AMES X and Ac R, MA is the set 
joo: xX € M and t e A). 
The phase map determines two other maps when either x or 
© is fixed. For fixed t in R the map 1 :X > X defined by 
a) called a transition; for fixed x ¿in X the 


map T:R > X defined by T (t) =t is Called a motion 


through x. 


MOREM 2.5 For each t in R, n’ is a homeomorphism on X onto 


itself. 


BRGOOE For any t in R, ne is clearly continuous because T 
tS continuous. m 2swinleetive since xt = yb implles x= y 
because x = x0 = x(t - t) = xt(-t) = yt(-t) = y[(t - t) = yO 


=y. To see that m fmoumyectave, let y € X be arbitrary. 
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í 


Then n” (x) = y where x = y(-t). Finally to show that ae 


MaS a continuous inverse, it is sufficient to show that 


ee or Ts lo see this we note that if 7° 


S a $ iy 
Sem are two transitions, the composition t on? 


CET S, 


is the 


transition T LOrE NYAK in 2, 


neon? (x) = m (n’(x)) = 1 (as) = xs(t) = x(stt) = q 


O 
It is also clear that N is the identity transition because 
Mor tany x in X, n(x) = x0 = x. Now since 


non’ = a’ 7 © = un the transition n” is the inverse of 


t 
T 


REMARK 2.6 For any x in X and [a,b] ER, the set xla,b] 


is compact and connected. 


PROOF The set {x} is compact and connected in X and [a,b] 
is compact and connected in R. Thus {x}x[a,b]-is compact 
and connected. Then since xla,b] is the continuous image 


of this set, xla,b] is compact and connected. Ve 


REMARK 2.7 Observe that the transitions nv, ieee Onc 


commutative group with the group operation being the compo- 


( 


Anon of transitions. 


PROOF The closure property is trivial for if t,s are in R, 


S t + s 
T 


foto iS in R and non = as we have already observed. 
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Associative: moon - SEN a usa) 


_ u LOS Vv 
= q or = (nm or )or 


nm (m(x)) = n°’ (xs) = (xs)t 


Commutative: ES 


m e) = yers) = (xt)s = nr (xt) 
Bern (x)) = mon’ (x) 
Identity: The identity is n° since mon? . qero 
E O+t 
= = iN = T OT 
Inverses: The inverse of e TE stimec 
Pa O0 
T oTt = T e. 


We now give a precise definition of what we have been 


calling trajectories. 


2.7 DEFINITION The maps y(x), v(x), and y (x) from X into 


X 


mare defined, for any x in X, 


Cw) = date te RI, 


Y (a) PS te RR, 


y (x) {xt:teR}. 
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For any x in X, the sets y(x), E, and y (x) are called 
mespectively the trajectory, the positive semi-trajectory 
pm@emehe Negative semi-trajectory. Note that for any x in X, 
y(x) = XR. 

Iwo trajectories Yı and Yo Mave a point in common, 
then the trajectories are identical. This fact is guaranteed 
Dyeune existence and uniqueness theorem. For if x e Yı 
and x € Yo> emeszrcheorem states that there exists a unique 
Scion passing through the point x. Hence Y, F Yo: This 
means that two different trajectories can never cross each 
Brenn This is an important fact and is very useful in 
dlagramming dynamical systens. 

Another example of a dynamical system which will prove 
to be useful is a dynamical system defined on a torus. 


Consider a differential system in the plane 
X = f(X) 


where f satisfies the hypotheses of the existence and unique- 
| i 
ness theorem. Set X = Ge and f = Cx and assume that the 
| 2 
MmMnetions fa and f, are periodic with period 1 in each of 


the variables x and y. Then 
IS e y +1), 1 = 1,2 


We take the square {(x,y): 0.< x <1, 0 < y < 1} and 


identify the sides x = 0 and x = 1 and the sides y = 0 and 
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EN. Then the points (0,0), (0,2), Groond Ll) are 
identified. Then we project the plane onto the torus by 
identifying any point (x,y) in the plane with a point (x,y) 
of the torus where X = x (mod 1) and y = y (mod 1). Then 
the trajectories in the plane projected onto the ee 


trajectories of a dynamical system on the torus. 


+ 


Er en Er >= 





Tt is clear from Example 2.3 that if we pick any point 
x on the trajectory, xt is also on the trajectory for any t 
in R. In Example 2.2 if we pick a point x on any trajectory 
and consider the positive semi-trajectory, then xt belongs 
to rx) for any t € R. The same results can be seen for 
the negative semi-trajectory. This observation leads us 


naturally into the following concept Ciena ance. 


9.8 DEFINITION A set MG X is called invariant whenever 
xt € M for all x e M and t eR. 


It is called Mos ve Ly invariant whenever this holds with 
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4 
replaced by R and is called negatively“ invariant Irit 
holds with R replaced by R`. 


From this definition we see that a set M& X is invari- 
ant if and only if it is both positively and negatively 
invariant for if M is invariant, we let x be an arbitrary 
BE ment oIM. Then xt e M for all t e R and in particular 
for t e R, xt e M and for t e R, xt e M. Thus M is both 
positively and negatively invariant. If M is both positively 
paemieratively invariant, then for any t e€ re xt GoM “and 
opa € R , xt e M. Thus xt e M for all t e Rand M 
is invariant. 

This concept leads us to some interesting and important 


properties associated with such sets. 


THEOREM 2.9 Let {M,} be a collection of positively invari- 


ame negavively invariant or invariant subsets of X, Then 


Í- e ë we ee ee ee ee 


eN A is positively invariant if and only if the set 
lees neratively invariant. M is invariant if and only 


wish is invariant. 


PROOF Let the sets M, be positively invariant. Let A =UM, 


i 
and B = AM, . For any x € A, x e M, for some i. Thus 
xt e M, for all t e R’ since M, is positively invariant. 


Mien t e A for all t e RT, Hence A is positively invariant. 


met x e B. Then X € M, for all i. Then since each M, is 


positively invariant, xt e M, for all te R’ and all 1. 
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Baus xt e B for all t e RT and Bers epost vely=invarrant. 

The proofs for the negatively invariant and invariant cases 
are entirely analogous. Let M be positively invariant. 

Let x € X-M and t e R and suppose xt £ M. Since -t € ar 
mot) e M. But xt(-t) = x0 = x which implies x s M and 

we have a contradiction. Thus xt e€ X-M and X-M is negatively 
invariant. The proofs of the converse and the third part 


are entirely similar. E 


THEOREM 2.10 Let MC X be positively invariant, negatively 
true. 


io icmeLosure ot tie cl M, and ats interior, Fint M, 


Mave the same property, 


2. Each of its components has the same property. 


ROOP., Consider the case of invariance. Let xe cl M 
and t € R. Then there is a sequence {x} in M such that 
Xn Zo ny ebhecernvarrance of M, xb esM for allt zR 
and a n. Then by continuity we have that xt > le eee AUS 
A el M and cl M is positively invariant. The proofs for 
the negatively invariant and positively invariant cases 
are entirely analogous. 

Assume M is positively invariant. Then X-M and 
hence cl (X-M) are negatively invariant. Consequently 


int M = X - cl (X-M) is positively invariant. The proof for 


M negatively invariant is entirely similar. 
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2. Let M be invariant. Now M = UM, where the M.'s 
pee Gisjoint components of M. Suppose M; is, not invarlant 
on some i. Then there exists x e M, amest e R such that 
me E M, . Now M; is a maximal connected subset of M and 
AROS] is a connected set. Thus if xt e M, for some j, 

M, U M, is a connected set and M, and M, are not disjoint. 
Bovis. xt Z M, herrea | which imoliies that xt ¢ M. This 
contradicts the invariance of M and hence each component of 
Mis invariant. For the converse let each component of M 
be invariant. Then since M qa and each component is 
w riant, it follows that M is invariant. The proofs of 
the other cases are similar. 

Suppose M is invariant, then so is X-M. Consequently 


Ci M and cl (X-M) are invariant. Hence bdy M = cl M f\cl (X-M) 


is invariant and int M = X - cl (X-M) is invariant. J / 


For the following theorem we need to define y(M), y*m), 
and y (M) where M C X. We define these as we would expect; 
that is y(M) = {Mt: te R}, y (M) = [Mt: t e R'} and 


Aim lariy for y (M). 


THEOREM 2.11 For any x e X, the sets y(x), y (x), and 

y (x) are, respectively, invariant, positively invariant 
endmesativelyzinwardant, A set McxX is invariant, posi- 
PEA invariant, 9or negatively invariant if and only if, 


4 = 
respectively, y(M) = M, y (M) = M, or y (M) =M. 
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PROOF Let xt. c y (x) for t. e R and let t.e R. Then 


Jl 
xt, (t,) = x(t, + t,) = xt 


2 


3 3 = E, + to. Since 


tjt, € R , we have tz € R and xt, Sence (0) 


JE 
Tor t 


is negatively invariant. The proof of the other cases is 
Prreiy analogous. For the second part, let M be negatively 
invariant. Now y (M) = {Mt: t € R }. In particular for 

t = 0, Mt = M and MC y (M). Suppose y (M) Æ M. Then there 
exists a t, € R such that Mt, % M. Then there exists 

x € M such that xt, £ M which contradicts the negative 
invariance of M. Thus y (M) C M and y (M) = M. Conversely, 
suppose y (M) = M. Then ÍMt: t e R} = Mand Mt C M for 


allteR. Thus M is negatively invariant. The proofs of 


the other cases are similar. 


The following remarks are useful in understanding this 


material and prove some useful results. 


BREMARKS 2.12 


l. The sets X and @ are invariant. 


OOF Let xe X. Then xt e X for all t e R since xt £ 2. 
Thus X is invariant. Since Ø = X-X, it follows that ß is 


invariant. 


A Set MEX is invariant, positively invariant, 
ermesatively invariant if and only if for each x e M, 


respectively, y(x) G M, Y'(x) EM, or y (x) CM. 
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PROOF Let M be invariant. Then by Theorem 2.11 y(M) = M. 
Maus (Mt: t e R} = M. Hence for any x e M, xt € M for all 
ETRE. But {xt: t e R} = y(x). Thus y(x)& M. Conversely, 
suppose y(x) © M for each x € M. Then (xt: te R} c M anå 
Aor all t € R and all x e M. Hence M is invariant. 


The proofs of the other cases are similar. 


Consider once again Example 2.1. For the points C 
n= tl, t2, ... we see that O = w = 0 and w = -kf sin 8 = 0. 
This shows that each x = Ce) does not change position with 
imrereasing or decreasing time. Such points are called 


= 


eati cal or rest points. We now state the formal definition. 


Poo rtON 2.13 A point x € X is said to be a critical 


point mix =) Xo for all t € R. 


REMARK 2.14 If x is critical and y # x, then x does not 


belong to the trajectory y(y). 


PROOF Suppose x is critical and y # x. Then {x} = y(x) 
pag@esince y(x) is invariant, {x} is invariant. Thus 
Pete) is invariant. Since y # x, y € X-1x) so y(y) € X-1x). 


Maus x ¢ yy). // 


Due remark, Figure 2.l should more properly be 


shown as follows: 
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nd 
SY > ae 


The next theorem contains several characterizations of 
fogrcal points but first the following lemma will be found 


useful. 


MENMA 2.15 If xe Xand x = xt for some t e R, then 


x = x(nt) for all integers n. 


PROOF If x = xt, then x(-t) = xt(-t) = x(t - t) = x so we 
need only consider positive integers. We complete the 
Meco by Induction. Clearly if x = xt for some Ut € EN 
ben x = Xt). Assume x = x(nt). Then 


At) see + nt) = x((n + Deo). // 


THEOREM 2.16 Let x e X. Then the following are equivalent. 
nic pent X 1s critical 


2. The singleton {x} is y(x) 


> 





Enessimelleton {x} is Y (x) 
Mier sin elagon {x} is yY (3) 
Demone keon tx} is xla bD] for some-a < b 


N n = UY) 





. There is a sequence {t !, t > 0, t, > O with 


xX = xt, for each n. 


PROOF The equivalence of statements 1. with 2., 3., and 4. 
are trivial. We next prove the equivalence of statemnts l. 
SmO. If x is critical then statement 6 follows trivially 
Since x = xt for all t € R. Conversely, assume there is 
meseguence {t }, t > 0, t. + 0 with x = xt. for each n 

n n n n 
Malet t e R. Ift = kt, aa A A SA O 
xt = x(kt_) = x by the lemma. Otherwise since t, > 0 for 
eeimtegers n, t/t. e R and is located between two integers 
anda: +). Thus k < t/t < k_ + 1 and 
n n n n n 

A = + - - : 
gee. et < Ck + 13t, Let 6 min UK, Dt, CG Ken 
Since en > 0, there exists an integer m > n such that 
En < 6/2. As above there exists integers K and = a 
wo hat kt <t < (k + 1)t.. Since this interval has 
m m om m 

length (K + Dt, A US 8/2, this interval must be 


contained inside the previous interval since 


max IK + 1)t, DA = kt} < 6/2. Thus we have 


SES tC (i F Te < k Flt. 
nn m m m m n n 


Cenne in this manner, we have constructed a sequence 


{kt} with the property that kt, + t. Now by the 
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continuity axiom, x(x 4) + xt and since x(k t) = x for each 
Oy che lemma, we have x = xt. Thus x is critical. 

To prove the equivalence of statements 1. and 5. we 
em emote that if x is critical, then 5. follows trivially. 
On the other hand, suppose x = xla,b] for some a < b and 
t e R. 1f t = kt, for some integer k and t 


l 1 
then xt = x(kt,) = x by the lemma. Otherwise, since [a,b] 


able 


is closed, there is a sequence {tt} DA a bbl vith t »a. 
Consider the sequence tt, - al. Clearly (t - a) > 0, 
E - a) > 0 and x = xa = x(-la) by the lemma. Since 

t e [a,b], x = xt. Combining these last two statements 
we have x = xa = x(-a) = xt, (-a) = x(t - a). Thus the 


hypotheses of statement 6. are satisfied and x is critical. // 


Erom this theorem we know that if we follow a point of 
any trajectory in a dynamical system and discover that that 
Bosnt as at rest for any period of time, then the point 


remains at rest indefinitely. 


Continuing to seek motivation from our examples, we 
again refer to Example 2.3. We see that if x is any point 
on the trajectory, x = X(X,>t) for some t e R but also 
x = X(X_ ot ween) tow ma = il, 22, 985, Transfering this to 
our notation of a dynamical system, we have that 
x = Xot = Xt + 2nt). As we would expect, any point with 
this property is called a periodic point. Formally, we 


State: 
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BERTINITION 2.17 A point x e X is said to be perio cake if 


there is a T ž 0 such that 
Beer: + T): for allte R. 


Amber T for which this holds is called a PERCIO TOLE: 
Mas 2 periodic point then both the motion and the 


ajectory are said to be periodic. 


The restriction that T # 0 is necessary since it is 
mere that xt = x(t + 0) for all t but the point x may not 
Bemperiodic. It is also noted that any critical point x 
is periodic since every T € Ris a period for x. The 
following theorems and remarks prove some of the expected 


properties of these points. 


THEOREM 2.18 If {x } is a sequence of periodic points with 


positive periods T, > 0 and x, * x, then x is critical. 


PROOF As in the proof of Theorem 2.16, given any t € R, 

there are integers k and k, + 1 such that 
eee tee tT = kT FT 
nn - n n nn n 

Since Ta 22 08 ko Th =. sahen x = Xara since Xn is periodic 


and xT = KT) by lemma 2.15. Combining these last 


two statements we have E xT = x. (KT) -> xt. Now 
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since % * X we Dave xxt. As. TIEN as arbitrary, 
xX = xt for all t eR. Hence x is critical. // 

micom the definition we can see that a point X e X is 
Begeodic if and only if there is a T # 0 with x = xT. For 
mt xX is periodic, we know there exists a T Y O0 such that 
mont + 1) for all t e R. But in particular, for t = 0, 
a have x = xI. On the other hand, if there exists a T # 0 
Such that x = xT and if t is any element in R, then 
EE t) = x(T +t) = x(t + T). Since t was arbitrary, 
we see that x is periodic. 

In Example 2.3 we noted that T = 2171 was a period of x 
and T = Un was a period of x. Clearly in this example, 
x has an infinite number of periods. We would expect that 
x would have a smallest non-negative period, which we will 
call its fundamental period. If there exists an infinite 
number of periods for a point x in a dynamical system whose 
fundamental period is 0, then as an immediate consequence of 
mee last theorem, we know that x is critical. This is true 
Since we can certainly form a sequence oe Or sche eds 
such that T + 0 and since the sequence tx, ) where x, = x 
tar all n has the property that Xx, + x, the hypotheses of 
the theorem are fulfilled. The following theorem proves 
the existence of a fundamental period and also proves why 


we need only consider positive periods. 


POR STE xX e X is periodic but not critical, then 


Biere is ar> os uch that T is the smallest positive period 


a 





Be Further if T' is any other period of x, then 


N 
T' = nT for some integer n. 


PROOF Let P = {t > 0: t is a period of x}. P # Ø since if 
ea period of x, T # 0, then x = xT. Thus 

x(-T) = xT(-T) = x(T-T) = x so -T is also a period of x. 
e either T or -T is positive, P A Ø. Now set T = inf P. 
Such a T exists since P ís bounded below by 0 and hence 
Man intimun. Clearly T > 0 so suppose T = 0. Then there 
is a sequence tt) IA E, + 0. Since x = xt For 
wwe x MUSt be a critical point by Theorem 2.19. Hence 

TY > 0. We see that T is also a period of x since there is 
AsSequence tt) a la > > T because T = inf P and 


x = xt oud en. Thus Dy the continúity axiom, 


k xt Tr Thus T is a period of x and by definition 

ts. the smallest positive period of x. Finally let t e R 
PA period of x. If t # nT for any integer n, then 

there exists an integer n such that nT < t < (n + 1)T. 

But since T is a period of x, nT is a period of x by 

lemma 2.15. Thus we neue x = xt = x(nT) which gives 

xt(-nT) = x(nT)(-nT) = x(nT-nT) = x. Then 

Mt (=nT Y = x(t-nT) which shows that t-nT is also a period 


of x. But since O < t-nT < T we have a contradiction 


since T was the smallest positive period of x. 7 


+ 
REMARK 2.20 If x € X is periodic, then y(x) and y (x) 


are compact. 
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| PROOF Let x be periodic with positive period T. Then 
Ax) = x[0,T] and v(x) IE TI cierto any t e LO. TI, 
EE OT] and for any t with t > T, t = t, + kT for some 
integer k and t, € TICs similar results hold fer 
momo. From Remark 2.6, x[0,T] is compaet. Thus y(x) and 


(a) are compact. // 


mom this analysis it is now clear that if any trajectory 


A A mm Se eee Tee Se 


trajectory must be periodic. For suppose xt, = xt, for 


Er < ta. 


weens shows that x must be periodic. 


Then = x(t, — ty) = dl ei Ie = t, = ty >02 
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MERO EONGATIONS AND LIMIT SETS 


erring again to Example 2.2, it was noted that in the 
damped spring case the solution converges to the point a 
as t + œ. This solution converges in the sense that there 
is a sequence ict with E, e R for all n, E, + o, and 
0 


xt, > (9): tee point is called a positive limit point. 


Formally, the definition is stated: 


DEFINITION 3.1 Define maps A", A from X into 27 by setting 
mr each X £ X, 

l. the set A (x) = {y e X: there is a sequence 
tt} MR Walt i te > e and xt, > wel: 

2. the set A (x) = ly e X: there is a sequence 


a m R with E, > -«w and xt > y). 


Morrany x e X, the set A (x) isc ica Sas ee 


set, and the set A (x) is called its negative limit set. 
REMARK 3.2 If x e X is periodic, then A (x) = A(x) = y(x). 


PROOF Let x £e X be periodic. By Remark 2.20, y(x) is 
memoact. Let ye y(x). Then y = xt for some t € R. Since 
mais periodic. y = x(t + T) = xT(t + T) = x(t + 2T). Clearly 
this idea can be extended inductively so consider the 

Beamence {geteml). Since T > 0, t + nT > +t and y = x(t + nT) 
werall n. Thus x(t + nT) > y and y € AT (x). Hence 


ae AT (x). Now let y € A’ (x). Then there exists a 
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sequence {tt} Ins with ol ace erallieny, 
xt, € y(x). Hence (xt } C y(x). Since this sequence 
Bonverges, it must converge to a point in y(x) because 
fo) is compact. Thus y e y(x) and Nixye y(x). Hence 
A* (x) = y(x). The proof for A (x) = y(x) is entirely 
analogous. 77 
An immediate consequence of this remark is that if a 
point x is a critical point, then A (x) = A(x) = x. From 
these two observations, it is an easy task to determine the 
positive and negative limit sets of Examples 2.1, 2.2, and 
Ms. Decause in Example 2.1 we have critical points, 
periodic trajectories, and trajectories which clearly do 
Mor converge. In Example 2.2, we have either a periodic 
trajectory for the undamped case or a trajectory whose 
Positive limit set is the origin and negative limit set is 
empty. Finally in Example 2.3 we have a periodic orbit. 
However, a point may not be periodic and yet have more than 
se point in its positive or negative limit set. The 


following examples thelp to illustrate this idea. 


EXAMPLE an Consider the differential system defined in 


Ro by the following equations in polar coordinates: 


ar p 

T~ r(l - r) 
dore 

mo 
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An extensive analysis is not necessary here because all the 
facts needed can be derived directly from the equations. 

This system does in fact define a dynamical system [Bhatia, 
Szego, p. 20]. Since the change in 8 is constant, only the 


change in r affects the system. Clearly for r = 0, 


a o Che origin is a critical point. For r= 1, 

= = 0 so r is again constant with time and the solution 
describes a periodic trajectory. Suppose r < 1. Then 
= is positive so r is increasing with time. If rol, 
then sg is negative and r is decreasing with time. such 


at 
trajectories are shown in the diagram below. It is noted 


that although — ISS positive Ior r < 1; P will never 
ar 


akh r= 1, for as r > l, a > Oh 





RIGURE 3. 2 


+ 
AS this diagram shows, for points p with 0 <r <1, A (p) 
is the unit circle and A (p) is the origin. For points p 


Ao e > l, A'(p) is the unit circle and A (p) = Ø. These 
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facts become clear when it is realized that a sequence 

{tt} can be found with > > © such that xt approaches any 
Pven point on the unit circle. Since the unit circle is 
Mode and the origin is a rest point, their positive 

ana negative limit sets are obvious. The direction of 

the trajectories is counter-clockwise because 8 is increasing 


with time. 


EXAMPLE 3.3 Consider the differential system defined in 


R< Berne diferential equations 


Era} 


>d o 
tt 


g(x,y) 


ó o 
N 


where the functions f and g are defined by 


ea 


Do) = E 
A E E 
AE YA OA 
O a 
g(x,y) =<4 -l if x < -1 


e ee 


Here the functions p(x) and q(x) are any continuously 


differentiable functions satisfying the following conditions: 
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ORCID) > a o A 
p(x) = 0 if x < 0, 
E yl 
DE x = 5 ES 
DE < l Dr 0 
ay 2 I en 3 
qa(y) = 0 O 


This system defines a dynamical system [Bhatia, Szego p. 21] 
and it is analyzed in the same manner as the preceding 
example. For x = y = 0, f(x,y) = 0 and g(x,y) = O so x and 
y are unchanging with time. Consequently, the point (5) is 
ical point. For x > 1, f(x,y) = 0 and g(x,y) = 1. 

ine trajectories of such points are straight lines parallel 
memume y-axis with direction in the ty direction. A similar 
ameltysis Shows that for x < -l, the trajectories are stralght 
lines parallel to the y-axis with direction in the -y 
direction. Now suppose 0 < x <1. Then g(x,y) > 0 so y 

is increasing with time. Bey < O.. f(x,y) > 0 sox is 
increasing. If y > 0, then f(x,y) < 0 so x is decreasing. 

As imitar analysis holds for -1 < x < 0. Thus the solutions 
fae cu be spiralling trajectories. .The question is do they 
Spiral outward or toward the origin or are they periodic 
trajectories. Such an analysis from the equations is very 
Mavo llved but the trajectories do in fact spiral outward 
[Bhatia, Szego, p. 22]. The trajectories are sketched in 


the diagram below. 
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URSS 


Brom the diagram one can see that for any point except the 
ieee in the strip -1 <x < 1, its positive limit set is 
the two straight lines x = +1 and its negative limit set 
Mene Origin. Since the origin is critical, its positive 
limit set is itself and for any other point in R°, its 


positive and negative limit sets are empty. 


Proceeding now with the concept of limit sets, some very 


important theorems are useful. 


THEOREM 3.4 For any X € X, 
l. the sets A* (x) and.A (x) are closed and invariant, 
2. the closure cl q) = o) y AT (x) and the 


altre cl (y (x)) = y (x) VU A (x). 


PROOF 1. Consider the case of A'(x). Let ty} be a 


sequence in A* (x) with U: If it can be shown that 
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+ + 
y € A (x), then A (x) is closed. Proceeding in this manner, 
it is noted that for each positive integer k, there is a 


sequence ee in R with ES > +% and xt + y, since each 


k 
Y, € A* (x). We may assume without loss of generality that 
ey, xt‘) < 1/k and E > k for n > k where d is the metric 
for the space X. These assumptions may be made since the 
xt“ get arbitrarily close to Yy and hence any point in the 
sequence at distance greater than 1/k from y, may be removed 
with the resultant subsequence converging to Vane The same 
reasoning applies to o A DA SACO Che Sequence 

{t }in R with t, = bo. The significance of this "diagonal" 


sequence is in consideration of {xt} and for purposes of 


clarity, we sketch the latter sequence below:. 


Yx Yo: A 


0 ==-—-= ve rn ee Y ers Y rs Y AT um m AA J 
* . . e e = * e* 
“ . SO a 


. t.’ 
G 
e@ @y 


s 
s’ oe 


foo) eee rss 
Pi LY YR ot 
A tt 
ae oe jf 
i: iad. ) 


© ® 
ger Keen 
Then Es Zend che claim is that xt, iyis TO See risa 


observe that 
d(y,xt,) E a(y,y,) + aly xta) < a(y,y,) + 1/n. 


Since 1/n and d(y,y,) tend to zero, this implies that 
+ 
aly,xt, ) > 0. Consequently, xt, > y and y e A (mm. Thus 


AT (x) is closed. To see that AT (x) is invariant, let 
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yE A’ (x) and t e R be arbitrary. There is a sequence 
im Mc E, + to and xt > Y. ¿men by theseentinulity 
axiom, xt (t) + yt. Since xt, (t) = x(t + t) and 
(t +t) > +» this implies that yt e A` (x). Thus A (x) is 
invariant. 

To prove the second part, it is recalled that MES = xR* 
By the definition of A (x) it follows that 
cl (y*(x)) D y* (x) U AT (x). Now let y e cl (y'(x)). Then 
there is a sequence Ur in MES such that Ya a: Since 
this sequence is in ES it follows that Ye > ee for some 
am € R. Hither the sequence it} has the property that 
meee tO, in which case y eh (x), or by the Bolzano-Weierstrass 


n 


theorem there is a subsequence (ti ) with ize + te R”. 
k k 


But then xt, > xt which is in ~ (x); and since also 
k 


xt > y, it follows that y = xt e bo. Thus 
K 


el (y*(x))C y (x) U A (x). The proofs for the negative 


case are entirely similar. // 


The properties of limit sets Just proved will be found 


to be very useful in further work. 
+ + + 
THEOREM 3.5 For x,y € X, if y e y (x), then A (y) = A (x). 


PROOF Let y € vy (x). Then there is at e R such that 
y= xt. Now let ue AT (y). Then there exists a sequence 


At Jin R with et such that yt. > u. Consider the 
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sequence {t + e) m ro early test Wee gen 
yt. > u implies yt. = xt (t_) = x(t + t) > uü. THis 
u € AT (x) and hence Mine O 

On the other hand, let ze N Can Then there is a 
sequence Kr) in Rwith t, > t+= and xt = z. Consider the 
sequence tt eu) in) Raeeclearly (t - t) > © since bt, > 
Then xt, > z implies xt, = x(t + E - t) = xt(t - t) = 


Wta - t) > z. Thus z € A’ (y) and N(x) ce A (y). 


ORK 3,6 Let x,y e X such that x is not periodic. Then 


the set Mty" (y): E y (x)) A 


BEROOF Since y e y (x), there is a t e R’ such that Ee ace’ 
Now suppose s e Mv (y): y € ES Then s € y (y) which 
implies that s = yt, for some E, € ay But since y = KE, 
it follows that s = xt (t,) = x(t + t,). Let t, e R such 


that t, > t + t,. Then for y = xt,, s £ uo Thus 

S £ Niy (y): y € y"(x)} and we have a contradiction. Hence 
+ ~ 

Wee y): y ey (x)=. Vi, 


REMARK 3.7 For any x F AT (x) =f\ {cl (y*(y)): ee vy (x)} = 


Nty' (zn): n is an integer). 


BROoOr If x is periodic, the proof is trivial since 

cl (y*(y)) = A (x) for all y. So we assume x is not periodic. 
Let {t} be any sequence in R such that t, + +». Without 

loss of generality it may be assumed that ba > 0 fog all n 

since any points in the sequence for which this does not 


hold may be removed from the sequence without affecting 
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ta > +, Since the intersection in the statement is taken 


, + 
E aell possible y's in y (x), set i xt, where each 


te {tt}. Now ge) = ty t: t e R} 
The claim is that N{cl (y"(y})} = N tel (y*(y,))}. 

By Theorem 3.4, fV{cl (vy (y))} = Nfy*(y) U A*(y)) = 
Niy’(y)}U Nat) ana 
Nta 0) + Nut) Va NN Na yy, 
By Remark 3.6, Ni) = f= Nix (y,)} and by Theorem 3.5, 
A" (y) = A*(x) = A (y_).  Thus 
Nte (vFty))} = Ni) = at(x) ana 


A tel Ce = Ni (y 0) = AT (x). Then, as claimed, 


{x(t tp E RT}, 


N {ei (y*(y))} = A {ci AI) and as shown above, 

Mier Y): Vane CO = AT (x). To prove the second 

part, let the sequence {it} = {In} Shen the proet Is 

identical. VL 
In the following proof and elsewhere in this thesis, 

the notation S(M,e) and S[M,e] will be used for the sets 

ala, M) < e, and {x: d(x,M) < e) respectively where 

Tre € Rand d 1s the metric for the space X. In the 

case where M = {x}, the notation will be shortened to 


meee and Sol x,c |. 


THEOREM 3.8 Let x e X. 1f£ cl (y (x)) is compact, then 


N (x) is a non-empty compact and connected set. 


PROOF Let el (y*(x)) be compact and let (t_) be any 


Sequence in RT such that un + +o, Now consider the sequence 
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{xt}. Since xt, e cl (yt (x)) for alin andee] (y* (x)) 
is compact, this sequence has a convergent subsequence 
{xt, }which Conver acewrema point y e el (y* (x)). Thus since 
ty + +o and xt, + y, we have AT (x) A Ø since y e A" (x). 
Now cl (y (x)) = y (x) U A’ (x). By Theorem 3.4, A’ (x) is 
eosed set and since it is a closed subset of a compact 
set, A (x) is itself compact. 

The proof of connectedness is by contradiction. Suppose 
A (x) is compact but not connected. Then O = PU Q 
where P,Q are non-empty, closed disjoint sets. Since 
cl (v’(x)) iezeonpaer, eis locally compact and hence there 
eistes an € > OÖ such that S[P,e] and S[Q,e] are compact and 
Seeyoint. To see that S[P,e] is ‘compact, it is noted that 
in a locally compact space, every point has at least one 
compact neighborhood. Taking the union of all these compact 
neighborhoods, yields a compact neighborhood of the set P. 
Then there is an e > 0 such that S[P,eJ is a closed subset 
of this neighborhood and hence S[P,e] is compact. The 
same reasoning applies to S[Q,eJ]. Since a metric space 
jis normal, there exist disjoint open sets containing P and 
Q respectively. Hence for a proper choice of e, the sets 
Ame] and S[Q,ce¢] are disjoint and compact. 

Now let y e P and z e Q. Then there exist sequences 
it} and CT with E, + +o and ae +o such that ae y 
and XT, > Zo. Without loss of generality it may be assumed 


that xt, € ses), XT € SIQ,Ee) and TI aaa ele 
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Since the trajectory segments xlt_,1,], A A a 
compact connected sets, they clearly intersect S(P,e) and 
S(Q,e). Thus in particular, there is a sequence ad 

with t, < Ea < T, Such that xT € H(P,e) where 

H(P,e) = {x: d(x,P) = e}. Since H(P,e) is a closed subset 
E el, H(P,e) is compact. Since ae is a sequence 
in H(P,e), it has a convergent subsequence {xT} which 
converges to a point y' and because Th > to, it follows 
nat y' € A (x). But y' ¢ PU Q which contradicts the 


EStslmption that a) = PU Q. Hence A* (x) is connected. // 


Since only the local compactness of cl y) was used 
in the proof, an immediate corollary to this theorem is 
Marat Che space X is locally compact, then A” (x) 1S Come eted 


Nr mever it is compact. 


The converse of the previous theorem holds if the space 
mers locally compact. However, the proof relies on the one- 
Pomme cOompactification of X. The fact that the one-point 
Gompactification is an extension of a dynamical system to 


enouhner dynamical system is proved in the following lemma. 


LEMMA Given a dynamical system (X,R,7) on a locally compact 
space X, let X' = X U {w} be the one-point compactification 


ete. and define T': X'xk > X' by 


Mt) TE EA: ER 
IAS & = 
W EA E nR 


61 





Men (X',R,7') is a dynamical system on X'. 
mmo First, by definition of T', 


nero) = lx ex 
m'(x,0) = 
W if x =o 


mags n'(x,0) = x for all x e€ X'. 


Beeendly,’if x € X, then m*(m*(x,t,380) = T'(x,ty + to) 


follows trivially so suppose x w . Then 


NUS 


T me ™'(w,t, + ty). 


+ ty) for all x € X'. It remains 


m'(m'(w,t,),t5) 


m'(m'(x,t,),t,) e 


1 
to show that m' is continuous. To see this, we will utilize 
the characterization that if the inverse image of an open 
set is open, then the map is continuous. 

eet U be any open set in X'. If w g£ U, then U is 
open in X and thus (m7) = (u). Now 72 (y) is open 
because T is continuous and hence (m't) TLU) iSv Open Ar 
we U, then X-U is compact and since w £ X'-U, X'-U = X-U 
and ce (nt )7+(x'~-u) - (nt) ~+(x-v) = UA Now since 
T is continuous, 1+ (X-U) is closed and hence Cr mex =p 
is closed so (mtu) is open. Thus the inverse image of 


every open set is open and hence T' is continuous. // 


THEOREM 3.9 If the space X is locally compact, then for 


=> me i A A A A A A AAA SS — Ru _— ë< 


any X e X, cel(y’(x)) is compact whenever N“ (x) is non-empty 


and compact. 
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PROOF Let A" (x) be non-empty and compact. Consider the 
extension of the dynamical system to the one-point compacti- 
@ecarion X' = X U {w}. Certainly, for any x e X, 

Eo) = A* (x) U {u}or ES - Ma) where ES is the 


positive limit set of x in X'. However since AT (x) is 


empact, it follows that Nr) Ne) holds because 


AT (x). Also we have that 


N* (x)=cly,(A*(x)) = 0, (A*(x)) 
r any x £ X, vee) = Y (x) since) = Ti) 


we have 


+ + + + + + 
Cl Er (9) = y (x) UN (x) = y (x)U A (x) = elyly (x)) 
and since AO x) is a closed subset of a compact space, 


fees compact. Thus cl(y’(x)) is compact. // 


THEOREM 3.10 If cl(y' (x)) Is conpaeer, then alxt,A’(x)) => 00 


as t > +», 


PROOF Let el(y’(x)) be compact and suppose d(xt ,A*(x)) AO 

as t gs“ Then there exists a sequence lad in R such CHEC 
ae te and d(xt 100) > eliñor all nio? N where N islan 
integer and € > 0. But since cl (y’(x)) is compact, the 
sequence {xt} has a convergent subsequence (xt, 3} which 
converges to a point p e cl Cy" (x)). Now since E, + +0, 

ae TOlLilows that E, SIA Since xt? p, we have that 

oe AY (x). Also since xt, > p, this sequence gets arbitrarily 


close to p and hence d(xt, ,p) < € for k sufficiently large. 
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| 


4 
Hence we have d(xt,, A (x)) < © for k sufficiently large 


but this is a contradiction and the theorem is proved. Var 


The study of limit sets can be very useful in predicting 


the behavior of a system. If the system describes a physical 


problem and the trajectories are sketched, it is a relatively 
simple task to discover the limit sets and hence predict 

the future behavior of the system. For example, consider 
again the pendulum problem, Example 2.1. If it is known 
that the pendulum is started at the equilibrium position and 
given an initial angular velocity of w = 2k, then one would 
Beedict that the pendulum would approach the position where 
6 = m and would never change direction to swing back. This 
result would be known without having someone watch the 
meauwlum to actually see what happens. It is: true that 

this is a very elementary problem, but the ideas involved 
have broader applications to problems of much greater 


Gatiiculty. 


It is important that the reader have a firm understanding 
Pace temiliarity with the symbols and sets thus far defined 
because the following additional sets and symbols become 


Aute important in the next two sections of this thesis. 


BERINITION 3,11 For each x £ X, define the following 


subsets D*, D, ae Io X, 


l. the set D* (x) = {y e X: there is a sequence eo 


in X and a sequence tt} in R such that x >» x and xt > as 
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2. the set D (x) = {y e X: there is a sequence {x} 
in X and a sequence {t} in R such that pe enim tS yr, 
3. the set Fa = ly €e X: there is a sequence ae 
in X and a sequence tt Jin R” such that ten 
and x,t, > ime 
4, the set J (x) = ly e X: there is a sequence {x} 
in X and a sequence (tt }in R such that ae, 
macy xt + yl. 
For any x € X, the set D’(x) is the first positive prolongation 
of x and the set D (x) is the first negative prolongation 
One X. 
The sets I” (x) and J (x) are called, respectively, the 


fans & bectutve and che first negative prolongational limit 


set of x. 


These sets are called "first" because higher-order 
Broloneatjons can be defined, However, in this thesis, 
such higher-order prolongations will not be discussed. 
Consequently, the word "first" will be suppressed. 

Rom, che definitions it is clear that any x e X, 
MES Y CO. DIO SY A) IRA (x), and I(x) DA (x). 
To achieve an intuitive feeling for these sets and to see 
Bnarsthe inclusions just noted may be proper, we present 


the following example. 


EXAMPLE 3.12 Consider a predator-prey problem where it 
is assumed that a single predator population feeds on a 


single prey population. Let x denote the size of the prey 
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population and let y denote the size of the predator 
population. Then this problem can be expressed by the 


system 


ax 
at — ax - bxy, 220, 
I = cxy - dy, y 2 0, 


where a,b,c, and d are all non-negative constants. 

These equations support our intuition since we would 
expect that any change in the size of the predator population 
depends not only on the size of the prey population but also 
on the size of the predator population. For example if the 
prey population is small and the predator population is 
large, we would expect that the predator population would 
Becrease, This expectation is confirmed by the equations. 

Reiner than solving this. system imstermssof t, it will 


be analyzed directly from the differential equations. 


Clearly for x= 0= y, = = 0 = SE and the point (0) is 


e ici cal point. One can also see from the equations 
2 ee erieleal polze FOr anye Pant 


= x 7 ax dy = s 1 
p (9); we have SE > 0 and at 0. Thus its associated 


trajectory lies on the x-axis. Physically, this is the 


that the point 


case when there is no predator population so the prey 


population increases. For any point p = Gok = 0 and 


= < 0 so its trajectory lies on the y-axis with direction 
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in the negative y-direction. This is the case when there is 


Mp rey population so the size of the predator population 


@ecreases. For small x and large y, = < 0 ana qe 0 
Thus both x and y are decreasing and x continues to do so 
WAtil y = a/b. At this point “= = 0mand as- y continues to 
decrease, x begins to increase. When x reaches a value of 
a/c, Sz = 0 and as x contínues to increase, y begins to 


increase. The physical situation described by this case 
maurice clear. The trajectories for this system are 


sketched below. 


Y 
do 
— _. 
mm — 
a y dle E 


BIGURE 3.172 


For any point p = (2) we have A (p) = de but 
I* (p) = Br: AO The fact that A” (p) = (0) means that 
if there is no prey population, the predator population will 
eventually decrease to nothing. The fact that Jt (p) = C) 
means that even though the prey population may decrease 
considerably so that the predator population decreases to 


nothing, there will be enough of the prey population left 
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to eventually increase its size to any amount. We also 


0 


Mebilce from this diagram that for any point p = os 


A* (p) 
D"(p) = y*(p) UH), D’(p) = y”(p) and J”(p) = 8. For 
any point p = (9), D"(p) = y’(p), I*(p) = 8, D’(p) = y (pJU 


Kon, AT p) = U}. For any other point p, 


(5) but J*(p) = ((%):x > 0}. Also 


- + ~ - + - 
me) = Y (PY, D (p) = y (p), J (p) = B = J (p). It may 
be argued here that we cannot form infinite sequences since 
"size" is not a continuous function. However, in an actual 
experiment, the following values for the constants a,b,c,d 
femewovoserved: aw~tl, b-~ .01,C - 3 x 1078 A 
These values were obtained by observing paramecia feeding 
= changes sign 
at 8 8 
is then at y ~ 100 and the value of x for which =- changes 
wn s at x~ 1.5 x 10°, These values are large enough 


Su east cells. The value of y for which 


to assume continuity. Also any mathematical model idealizes 
the physical situation: the question is how "close" [Bailey, 
Br 245, 246], 

New we consider some basic properties and important 


theorems about these limit sets. 


THEOREM 3.13 For any x e X, 
INES et D’ (x) is closed and positively invariant, 
Zo sen J*(x) is CloOscdmand inyanlant, 
3. the equality D(x) = 006 sx) Poldo: 


Analogous results hold for D (x) and J (x). 
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PROOF Statement 1 will follow from statements 2 and 3 so 


statement ¢ is proved first. This proof parallels the proof 


of Theorem 3.4. Let G be a sequence in I'(x) with 


Y + y. For each integer k, there are sequences{x\} Ir 


n 

k e k k k, k 
andit, }in R with x, * Xs, t, > te, and xt > y,. Without 
loss of generality, it may be assumed that t` > k, 


k | k, k 
< 
d(x ,x) < 1/k, and alx ata 


ame ee Ore Now consider 


the sequences {x}, ae Clearly ee E 


SS NW DS 


Fa too merc 


xt, weve This is true since 
Den 
may) < d(x toy) + aly,.y) < 1/n + daly y): Thus 


y € J” (x) and I” (x) is closed. 
To see that ES is invariant, let y € Be and tC ETR., 
There is a sequence {x} in X and a sequence {tt } in RÝ 


such that E, AS AAA xt. > y. Now consider the 


n 


Pe. 


sequence ho abi Cleary t, +t > +”, and 
x(t +t) = ate) > yt. Since x, + x, it follows that 
fe € a (x). Spice E R was arbitrary, it follows that 
J*(x) is invariant. 

Consider statement 3. Observe that 
E) Y(x)U Tx) always holds from the en So 
let y € D’(x). Then there is a sequence {x} in X and a 
sequence {t_} in EN with x Xs ve > y. It may be assumed 
that either t >te R* or t, >+, if necessary by taking 
subsequences. In the first case xt Fxtabythe continulty 


+ 
mom. Hence, xt = y € y(x) as t € R. In the secona case 
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y € Mix) Besefinition. Thus ye v(x) U J*(x) and 

Oc y (x) U ON Therefore, DR) = y (x) U J (x). 
Mo see statement 1, recall that cl (y* (x)) = y' (x) U eon 

But it has been shown that M(x) Ce J ee Hence 

D' (x) = cl (y*(x)) U I(x), because Ce cor Since 

Moon these sets are closed and positively invariant, it 


+ 
follows that D (x) is closed and positively invariant. Veh: 


The following theorem will also be found useful in 


subsequent sections. 


THEOREM 3.14 Let X be locally compact. Then A (x) # Ø 


+ 
whenever J (x) is non-empty and compact. 


PROOF Suppose A` (x) = Ø. Then 

el (y*(x)) = y™(x) U A(x) = y (x) since A (x) = Ø. Hence 
v(x) is closed. It is claimed that v(x) ds aisjoint From 
Fi) mor if yi (x) N J+ (x) # 8, then by the invariance of 
TO, Y'(x)c ECK Since Tx) is compact, for any 
sequence tt, in RË such that t, > to, the sequence {xt} 
must have a convergent subsequence. It follows that the 
point to which that subsequence converges is a point in 

A* (x) and hence AT (x) Z Ø. Thus Gon Ft = ğ as 
claimed. Since SOS is non-empty and compact, 

ae 0.7) = 0. -Thus since- -X i3 locally compact, there 
is a ô > 0 sueh that Sea is compact and disjoint 
from x), Now choose any y € Fi. There is a sequence 


+ 
{x} in X and a sequence {t,} in R such that x, > x, 
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| Es Dre, and xt + Y. it may be assumed that 
| X g Sid" (x),6], xt esti) ,6] for allen Then the tao- 

| jectory segments x, L0,t,] intersect H(I"(x),6) where 

| it’ (x),6) = {ze X:d(J* (x) ,z) = 6}. Therefore there is 

a sequence mn 0 < aS t> such that Sl HCI” Come Je 
Since H(J* (x) ,6) is compact, it may be assumed that 

ma, 7 2 CE (x) 5). By taking subsequences if necessary, 


it may be assumed that either cars Lee RT ORT Foe. If 


Ii 
N 


> tre Fi wen Dy the continuity axiom Kaas sa 
OZ € rx) which contradicts GON Slo (Coen = Y. 

If o to, then z € Con but thiszgeontzadTers 

Z € MO (x),5) as Jt) N I (x.s) = Y. Thus the original 
assumption that A* (x) = Ø is untenable and hence 


A* (x) #8. Tf 


THEOREM 3.15 Let X be locally compact. Then Fi 1 


Den empty and compact it and only if D* (x) is compact. 


PROOF Let IT (x) be non-empty and compact. Then by the 
preceding theorem, AT (x) is non-empty and, since it is a 
closed subset of cor eis ceampaet, FButztchenzgei (y’(x)) 
is compact by Theorem 3.9. Hence 

D’(x) = (0 U a*(x) = cl (y*(x)) U J (x) so that D(x) 
is compact. On the other hand, suppose D’ (x) is compact 
and assume ee ERO R But J" (x) = Ø implies vy" (x) is 
compact. Then y(x) = cl (y’(x)) so cl (y*(x)) is compac 


but this implies AT (x) # Ø by Theorem 3.8 and the assumption 


{1 





+ 
that J (x) = Ø is contradicted. Thus J*(x) # Ø and since 
+ 
J (x) is a closed subset of D’(x), it follows that J”(x) 


as compact. Y) 


Before proceeding to the next theorem, the set D’ (M) 
must be defined where Me Xx. Thus D*(M) = {y e X: there is 
a sequence Ea in X and a sequence tt} in R such that 


x, >» xeM, and nn > y). 


THEOREM 3.16 For any compact set MC X, the set D'(M) is 


closed and positively Twa ran ta 


PROOF Let ee be any sequence in D? (M) such that Y ye 

Then for each Y% there exist sequences (x/) in X and (69 

in Re such that xf > x and ir + y, Where x. ce M. MWEtnout 
n K nn k k 


loss of generality it may be assumed that dt y) < 1/10 


Now consider the sequences ag in M and (6) in R’. Then 
ZA Y, and since M is compact, it may be assumed that 
=: E TOUS pr AS AN Ce 
n nn 
nn nin 
d(x toy) < Rey) + OS) ee, d(y,»y)- Then 
by definition, y e D'(M) and hence D'(M) is closed. To 
< 
see that D (M) is positively invariant let y e D (M) and 
let t € EN be arbitrary. Now y € D’ (M) implies there exist 
+ 
sequences ae ee ond ee in R such that x »xeM 
+ 

and 2 > y. Now consider the sequence us +o) on he 

= + . 
Since ce > y, it follows that x_t,(t) US t) + yt 
Inús yt € D? (M). Since t e R was arbitrary it follows that 


D? (M) Ie tant. // 


fe 








We conclude this section with an important theoren. 


THEOREM 3.17 Let x,y in X. Then y e J'(x) if and only 
mex ce J (y). 


PROOF Let y € lx). Then there is a sequence lx) ire XK 


and a sequence {tt} in R such that x, = x, t, + +#, and 


u Set xt = yo» and EN on Then Yan > Y» 


Ba >, and ES Ne) = xt, - e er 
Consequently x e J (y). The converse holds in the same 


manner. Vy 
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Ove RECURSIVE AND DISPERSTUE CONCERES 


Using the sets defined in the last section, we continue 
with our analysis of dynamical systems by studying recursive 
and dispersive concepts. To motivate the idea of recursive- 
mess, consider again Example 3.3. Suppose we consider a 
set A to be the union of two circles centered respectively 
at the points 2) ana G with radius r where 0 < r < l. 

In this example, all trajectories in the strip -1 < x < 1 
spiral outward from the origin with the lines x = -l1 and 

x = l] as their positive limit set. Thus for any point in 
the strip -1 < x < 1, its trajectory will eventually 
intersect A. The trajectory will pass out of A again but 
it nevertheless will return some time later and repeat the 


meocess. This is shown in a following diagran. 





The set A is said to be positively recursive with 


respect to the trajectory. We now state the formal definition 


Of recursiveness. 
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DEFINITION 4.1 A set AC X is said to be positively 

meeursive with respect to a set BC X if for each T e R 

mere is a t > T and an xe B such that xt e A. Negative 
Mecursiveness may be defined by using the inequality t < T. 
Set A 1s self positively recursive whenever it is positively 


recursive with respect to itself. 


The simplest example of a self recursive set is the 
set {x} where x is a periodic point. By placing restrictions 
on the system by demanding that the entering of points from 
Binto the set A happen with some regularity, other concepts 
such as almost periodicity can be discussed although such 
is not our purpose here. We now define some concepts and 
establish properties in connection with the concept of 


recursiveness. 


DEFINITION 4.2 A point x e X is said to be positively 
Poisson stable if every neighborhood of x is positively 


recursive with respect to {x}. 


This definition means that the trajectory through 
the point x eventually intersects every neighborhood of x. 


ime following diagram illustrates such a situation. 
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To give a better understanding and some characteriza- 
mons Of positively Poisson stable points, the following 


theorem is presented. 


THEOREM 4.3 Let x e X. Then the following are equivalent. 

elie Pont iS positively Poisson Stable: 

E NEN a neighborhood U of x and a T >20, 
De DeNerxU eV) Tor some t >T, 

Br e pony xX is contained in lis positive “any 
set A'(x), 

4, The set cl (Y (0) = A'(x), 

5. The set y(x)E A'(x), 

6. For every e > O there is at > 1 such that 


ELSA. E 


PROOF The equivalence of statements l: ana 2 follows from 
ne definition and our previous discussion. To see the 
equivalence of statements 2 and 3, consider a sequence 

=} in Rt such that ue + +o, To each n associate the 
neighborhood S(x ,1/n) of x. Then for each t_ there 
exists at, >? m such that xt, € S(x Imn). Clearly 

> ero and xt, Jee. CnwSs xc AT (x). Conversely, let U be 
a neighborhood of x. Since x € AY (x) there is a sequence 
it} in R' with t > +> such that xt, » x. Thus xt, must 
eventually be in every neighborhood of x so that, in 
Barticular, given any T > 0 there exists a a > T such that 


xt, e U. Next, we establish the equivalence of statements 
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El, Let x € A* (x). Then since AT (x) is positively 
invariant, we have v(x) Nies Since 

cl (y*(x)) = y(x) U Mo, it follows that cl (y" (x)) = A(x). 
Conversely, let cl (y"(x)) = AT (x). But since we know that 
el (y’(x)) = Y'(x) U A* (x), it follows Chat (e A* (x) 
and hence that x € AT (x). To see the equivalence of 
Seavements 3 and 5, let x € A* (x). Then, since A* (x) is 
Invariant, xt € N for all te R: Thur ye AT (x). 
Conversely, let y(x) c A* (x). Then for any t e R; 

Rt e AT (x). In particular, for t = O we have x € AT (x). 
Finally, we prove the equivalence of statements 3 and 6. 

To this end let x € A*(x) and let € > 0. Then there exists 
a sequence {t } in RS with Le > © SWeh that xt, O 
Without loss of generality, we may assume E, O Ei A 
Then since xt + x, this sequence must eventually be in 
every neighborhood of x. Hence there exists a ta Ze 

such that xt, Br ec Conversely, choose a positive 
Del) sequence te} such that €, * 0. Then for each n, we 
have a t, 2 i such that xt, € S(x,e,) . Since €, >0, 

we have that xt, > x. Consider Che sequence {t 3 in EN 
Either en >to Or it contains a convergent subsedmence 


ft ) such that to >teRi. Ift >+ to, xe A(x) by 
N ny n 


definition. If t + t, then xt > XG. PUCE Ince E NR; 
ny ny n 


we have x = xt and therefore x is periodic with period t. 


ite totlows that x € AT (x) and the proof is complete. // 
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From this theorem it should be clear that whenever a 
point x is positively Poisson stable then so is xt for every 
me KR. This is Mts Dina if x IS positively Polssen 
stable, then x e A'(x) which implies xt e A (x) for all 

t e R by the invariance of A*(x). But A (Go) = A (xt) for 
Plt e R by Theorem 3.5. Hence xt e A (xt) and xt is 
positively Poisson stable. 

According to the above theorem, a point x is positively 
Potsson stable if x e A (x). SIN lA pointi IS 
defined to be negatively Poisson stable if x e A (x). It 
memsaid to be Poisson stable if it is both positively and 
mMesatively Poisson stable. If a point x is Poisson stable, 
then both its motion and its trajectory are said to be 
Poisson stable. 

From the preceding theorem observe that if 
cl (y'(x)) = A* (x), then x is positively Poisson stable. 

Now consider the situation where y* (x) = A’ (x). In chat 
se une point x 1s positively Poisson stable since 
cl (Y (x)) 


that (x) 


(x) U A(x) = A(x). We will now establish 


A" (x) same Only if x sea pemodle porn 
NOR x 15 a periodic point, it is clear that v(x) = A’ (x). 
On the other hand, if v(x) = A (x), then x € A* (x) and 
by the invariance of AT (x), y(x) = AT (x). Then 

KE! € Y'(x) RS OC and therefore therelexists aTr 


wen Chat xt’ = xf. Hence xt'(t) = xt(t') = xT(t). 
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Therefore xt = x(t + T - t') for all t e R, showing that 
w periodic with period T - t' > 0. 

With this last fact it is a natural question to ask 
whether there are Poisson stable points which are not 


periodic. The following example shows that such points do 


exist. 


EXAMPLE 4.3 Consider the dynamical system defined on a 


torus by the system 


at = Dee 


E 


qe A 


where f(x,y) is periodic in both x and y with period 1. 

ee is f(x,y) = f(x + ly +1) = f(x + lly) = f(x,y + 1). 
Berther, we have that f(x,y) > O for x and y not both 0 mod 1 
and f(0,0) = 0. The constant a is an irrational number. 


Clearly the point p = m is a critical point because 


ales point = = 0 = SL. Solving this system simultaneously 
we have = l/a which implies y = ax + c where ce is some 


wn ant offintegration. Thus the trajectories describe 


straight lines of slope a in the plane. We may assume 


a > 0. Now no trajectory can contain p except y(p) = {p} 
; . ax dy i 
since p is a rest point and since an and ar are both positive 


nor alk other points. 
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Anote now that the function y = ax + c is not defined 


Hache point p since at that point n = >. So we claim 
that there is only one trajectory which approaches the 
Point p as t +o. The following diagram will help us to 
determine this trajectory. For purposes of this diagram, 


we assume 0 < a < l. 


Y P 


Consider the trajectory y = ax + (1 - a). As x > l, It 

is clear that y * 1. Hence this trajectory approaches 

the point (2) as t > ©, But on the torus, the point 

p = (0) = Er. Hence this trajectory approaches p as t > œ. 
Generalizing this approach to the case of an arbitrary 
positive value for a, the trajectory which approaches p as 

t > » is given by the equation y = ax t (= (a mod ia 

To verify this fact, let x +1. Thenyrart (1 - (a mod 1)) 
=1+ (a - (a mod 1)). But a- (a mod Deals mar sone 
nteper n and hence as x > l, we have y > n + 1l. Thus 

as t > œ, this trajectory approaches the point E | D but 
Main, in the torus, this point and the point p are identi- 
fied. Thus this is precisely the trajectory which approaches 
meas > > ©. Cali this trajectory Y5: Similarly there ac 
eny one trajectory which approaches the point p as 


t is given precisely by y = ax +0.. It is clear 
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maat as xX > 0, we have y > 0 and this verifies that this 
mreajectory does in fact approach p as t + -~. Call this 
trajectory Yı° 

Now we claim that no trajectory in the torus is periodic. 
Tf there is such a trajectory, then it will pass through 
the same point twice and without loss of generality, we 
may assume that this point lies on the y-axis. We may 
make this assumption because the trajectory must intersect 
the y-axis. Assume that this point is the point ($) where 
we < 1. Then if the trajectory given by y = ax + c 


entersects this point twice, we have ($) = ( „mod a 


ae > 


an 
(an + c) mod l. 


meone n >00. But this implies that c 
Since 0 < c < 1, we know that c = c mod 1 = (m + c)mod 1 
where m is any integer. Hence it follows that 


k for 


(an + c)mod 1 = (m + c)mod 1. This implies that an 


some integer k. Now since n # 0, it follows that a k/n. 

But this implies that a is a rational number which is a 

eontradiction. Thus no trajectory can pass through the 

same point twice and hence no point except p is periodic. 
We also note here that the trajectories Ya and Y> are 

not the same for if they were, then the equations 

y = ax + (1 - (a mod 1)) and y = ax + O must be equal. 

But this implies that a mod 1 = 1 and thus a must be an 

integer which contradicts the fact that a is irrational: 


We now sketch several of the trajectories in the 


following diagram. 


Gall 











To complete our analysis of this example, we see that 
A (y) eipheand that u = {p}. For the point p, we 
N (p) 
A (y) 


+ 
have A (p) lp). For any other trajectory, we 


have My) 


the torus since these trajectories 
mot perlodic. Also, u = A (Yo) = the torus. 
Hence we have that all the points on Y, are positively 
Poisson stable since for each x € Ya XE Fu. Also 
points on Y> are negatively Poisson stable. All other 
points are Poisson stable. 

If we change this example so that f(x,y) > 0 Tortall 
mama y, then we will have no rest points and no periodic 
trajectories. Hence every trajectory is dense in the torus 
and the positive and negative limit sets of each point is 


the torus. 
We now introduce the notion of a non-wandering point. 


DEFINITION 4.4 A point x e X is said to be non-wandering 


every nelghboorhood U of x is self positively recursive. 
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Ms clear from this definition that if a point is 
PF positively Poisson stable then it is non-wandering. This 
idea can also'be seen in the following theorem which 


gives some characterizations of non-wandering points. 


THEOREM 4.5 For any x € X, the following are equivalent. 
i. Ihe point x is non-wandering, 
ANC no xls an clement OF LUSspoOsiiuve 


prolongational limit set ON 
peewee Mel ehnbomnood Of x 1s seli negatively 


reeursive, 


prolongational limit set J (x). 


PROOF Assume x is non-wandering. Consider a null sequence 
le.) 0 < E O, and a sequencett, ) o omc 
mce each S(x,s_) is self positively recursive, we have an 
n 

x, € S(x,€,) and at, > E, with x Th € S(x,E, ). Since 
IS 0 we have x > x and x. T. > x and since T_ > to we 
n 6 en n 
tonclude that x € iGo Thus statement 2 holds. Assume 
X € T(x): Then there is a sequence B in X and a 
sequence cae WB une 2 x and te such tha. 
xA > x. Now for any neighborhood U of x and T > 0 
a such that t > T, x. _ Ee U, and xt € U for 

i n n an 
mee N. Thus U is self positively recursive and, consequently, 
x is non-wandering. The equivalence of statements 3 and 4 


is proved in the same manner. It remains to show that 
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statements 2 and Y are equivalent. This follows directly 


irom Theorem 3.17. E 


We shall now prove some theorems which establish the 
connection between Poisson stable points and non-wandering 


Points. 
MABOREM 4.6 Let xeX. Every y € AT (x) is non-wandering. 


PROOF It is necessary to show that if y € A(x) for some 
+ 

EA, then y € J (y). Ifye A* (x), then there is a 

sequence {tt} do A o incer O 


we may assume, 1f necessary by taking a subsequence, that 


tatl = Un Or each n. Then setting 

MS Ea and xt. = x, we have x. + Y, x ta T ta) = 
xt 41 = vy, and GR tel ~ pa >» +0, Thusy €e J Gee and 

y is non-wandering. // 


THEOREM 4.7 Let PC X be such that every x e P is either 
positively Gunccocively forsson Stable. then every 


bee Cl P ES! non-wandering. 


PROOF Let {x} in P and o We must prove that 
XE Fi. For each n we have that either x € N“ (x) or 
x, € A (x_)- Thus by taking a subsequence if necessary, we 
may assume that x_ € AT (x ) forall nor that a Cae 

n n n n 
Malin In the first case, for each n there is a 
cm AAN d(x,»x,t,) < 1/n. Then clearly 


d(x,x,t,) < d(x,x,) + d(x ‚xt, ) 3 d(x,x,_) + l/n. This 
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shows that a Pand consequently z € Fu. In the 
E ond case similar considerations show that x e J (x). 


Thus by Theorem 4.5, every x € cl P is non-wandering. // 


A partial converse of this last theorem will now be 
presented without proof. The proof is entirely constructive 
Puee nov particularly instructive. A proof may be found 


in Bhatia and Szego. 


THEOREM 4.8 Let X be complete and let every x e X be 
Ber wandering. Then the set of Poisson stable points P 


is dense in X. 


To see how several of these concepts are related, we 
W Sent the following summary. It is clear that a critical 
wis a periodic point. Also, all periodic points are 
Poisson stable since the positive and negative limit sets 
of a periodic point are the periodic trajectory. Certainly 
miomperitodic point is contained in its periodic trajectory. 
Finally, all Pelee stable points are non-wandering points 
Boints for if x is a Poisson stable point, then 
x € ¡MES J* (x). These ideas are summarized in the 


following diagram. 









periodica P oints 


stable p oints 


Poisson 









Non-wander ing Points 
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Pemmew define a different type of stability. This 
notion will become important in discussing "dispersive" 


eoncepts. 


DEFINITION 4.9 For any x e X, the motion mn, is said to 
be positively Lagrange stable if cl (y’(x)) 1s compact. 
Further, if cl (y (x)) is compact, the motion ", is called 


Merpavively Lagrange stable. It is said to be Lagrange stable 
mel (y(x)) is compact. 


If X = R”, then the above statements are equivalent to 
the condition of the sets io, Y (x), y(x) being bounded, 
Meeti velly. The following statements about this notion 
are easily proved using Theorems 3.8 and 3.9 and Remark 3.10. 

Too Se lOcalily compacta thena motion La lus 
positively Lagrange stable if and only it 
AT (x) is a non-empty compact set. 

eTa moton ue 1S positively Laecerangene taole; 
then AT (x) is compact and connected. 

MOL Lon ur is positively Lagrange stable, 


then a(xt,A'(x)) + 0 as t > +o, 


We are now ready to discuss dynamical systems which are 
mMored by their lack of the properties of recursiveness. 
One can easily draw a parallel between these dispersive 
concepts and the recursive concepts just defined. First, 


we need some definitions. 
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DEFINITION 4.10 Let x e X. 


ur 


The motion ul is said to be Positively Lagrange 


unstable whenever cl (y*(x)) is not compact. 

The motion is said to be negatively Lagrange 
unstable whenever cl(y (x)) is not compact. 

The motion is said to be Lagrange unstable 

if it is both positively and negatively Lagrange 
unstable. 

itespointg x as called positively Poisson unstable 
whenever x £ A(x). 

ISO es called negatively Polsson tasa Te 
whenever x £ A (x). 

The point x is called Poisson unstable whenever 
IIS pothposSitively and negatively Eorsson 

Wine wie. 

Demont x 15 4called wandering whenever 

x £ Fi. 


Mats 1s a racher lengthy definition but the definitions 


are what we would expect. We are now ready to define 


dispersive concepts on the whole space X. 


DEFINITION 4.11 The dynamical system (X,R,7) is said to be 


re 


Lagrange unstable if for each x € X, the motion 
ae is Lagrange unstable, 
Poisson unstable if each x em is Poisson unstable, 


completely unstable if every x € X is wandering, 
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| 4, dispersive if for every pair of Points Ve 
| 


not necessarily distinct, there exist neighbor- 


hoods Ue ot x and U, of y cuch CRAT Sr ES AOL 


y 


positively recursive with respect to Uy and Uy 


is not positively recursive with respect to U. 


One may wonder why such concepts as positively and 
negatively wandering points are not defined. This is 
unnecessary since x e I) if and oniy ifr x euer), Ore 
can see from the definition that a completely unstable 
system must be Poisson unstable since if x £ Fo, then 
x £ Fo. The following remark proves the relationship 


between all the sets defined. 


REMARK 4.12 In definition 4.11, each statement implies 


the preceding statement. 


PROOF Suppose (X,R,1) is dispersive. Then for every pair 

Of points x,y e X, there exist neighborhoods U, of x and 

Uy om y such that UR is not positively recursive with respect 
to Uy: Set x= y. Then there exists a neighborhood of x 
which is not self positively recursive. Hence x is wandering 
so x ¢ Fu. Since this is true for all x “e X, We Paye 

that X is completely unstable. Next suppose X is completely 
unstable. Then x £ Fx) for all x £ X. Then also 

w (x) for all x € X. Hence x£ A*(x) and x £ N(x) 

since Mx) a Fi Dat ES Is Das ins Toison 


unstable. Finally assume X is Poisson unstable. Then 
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x £ A* (x) and x Z A (x) for all xe X. Suppose cl (y"(x)) 
is compact. Then by Theorem 3.8, AT (x) is non-empty and 
compact, hence complete. Then by Theorem 4.6 every point 
in At (x) is non-wandering and by Theorem 4.8, the set P of 
measson stable points is dense in AT (x). Thus Ps none 
empty. But if q is a Poisson stable point in Fo, then 
q € A*(q) Mene Space 1S not Polsson unstable.. Hence 
Nave a contradiction so that cl (y’(x)) is not compact. 
A similar argument shows that cl (y (x)) is not compact. 


Hence X is Lagrange unstable. // 


The following diagram provides a summary of the ideas 


in the preceding remarks. 







dispersive 







Camp letely unstable 


Poisson unstable 


Lagrange unstable 





Thus we know that a dispersive system is completely 
unstable, a completely unstable system is Poisson unstable 
and a Poisson unstable system ís Lagrange unstable. One 
may wonder if Che reverse implications hold. In general 


they do not as the following examples illustrate. 


To show that a Lagrange unstable system may not be 


Poisson unstable we want x e€ A (x) buvzal (y*(x)) not 
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compact. Since x € ON we know that y (OE AT (x) and 
hence cl a) = ED Thus we must have an example 
where AT (x) Pome umecompact and yet x € A’ (x). To this end, 
consider Example 4.3 of the dynamical system defined on the 
merus. We delete the point p from the torus so that the 
merus 1s no longer compact. Then for any x in the torus 
Bech that x É Y>, we have Ns = the torus and hence 
x € ¡MES but AT (x) is not eonpact. Mus cl (y*(x)) 
EOC compact. Clearly for x € Yq> A (x) = @ and hence 
ICY (x)) is not compact. Similar results hold for x € Yo 
and cl CUCID is not compact. Thus for all x in the torus, 
cl eo) is not compact and cl (y (x)) is not compact. 
Hence this system ís Lagrange unstable but it is not Poisson 
unstable because for any X £ Yo; xe A’(x). 

To see that a Poisson unstable system may not be 
completely unstable, consider the following phase portrait 


in the x,y-plane. 





inemunit circle contains a rest point p and a trajectory y 


+ — 
such that for each point q e y we have A (q) = A (q) = {p}. 
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mie trajectories in the interior of the unit circle have 
the same positive and negative limit sets as y. All tra- 
ieevories in the exterior of the unit circle spiral to 
Mmoemunacecircle as t + +>, so that for each point q in 
the exterior of the unit circle we have N (a) = {pI y; 
and A (q) =Ø. If we consider the dynamical system obtained 
won Chis one by deleting the rest point p then the resul- 
tant system is Poisson unstable. However it is not 
completely unstable because for each q e y, I*(q) = y 
Bo dq Ee no 

Finally to see that a completely unstable system may 
not be dispersive consider the following dynamical system 


defined by the differential equations 


ms 
ll 


Sin y 


2 
y = cos y 


This system contains trajecrorres Yk given by 
mer (X y): ve key. KEI Et 

These are lines parallel to the x-axis. Between any two 

consecutive Y] 5 the trajectories are given by 

y = {(x,y): x + ¢ = sec y), where c is a constant depending 


on the trajectory. The phase portrait between the lines 


y = -n/2 and y = 1/2 is shown in the diagram below. 
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Ya Vo 


can be seen that this system is completely unstable 
since for any p € Y_1> I*(p) = a ana forgall other points 
Benot belonging to any Yi, Jt (p) = %. Thus p £ J(p) for 
pay Pp and hence every point is wandering. However if 

pe Y_] and q € Vee then every neighborhood of q íis recur- 
Sive with respect to any neighborhood of p and the system 


is not dispersive. 


The following result shows that to prove Lagrange 
Mescebility, it is sufficient to prove that the space 


mompositively Lagrange unstable. 


THEOREM 4.13 A dynamical system (X,R,m) is Lagrange unstable 


Band only if it is positively Lagrange unstable. 


2200er The necessity is trivial. On the other hand, let 

the space be positively Lagrange unstable. Then for each 
mee X, cl (y*(x)) is not compact. Suppose cl (y (x)) is 
compact. Then A (x) is non-empty and compact. Since A (x) 
mewinvariant, yt € A (x) for allt e R? and ye NGO Tmn 
emee A(x) is closed, cl (y (y)) C A (x) and hence 


el (y* (y)) is compact. This is a contradiction. Hence 


Eje 








cl (y (x)) is not compact for each x e X. Thus the space 


is also negatively Lagrange unstable, yielding the result. // 


For a dynamical system defined in the x,y-plane, the 


following theorem provides an important fact. 


THEOREM 4.14 For dynamical systems defined by a system of 


BE rentilal equations in the euclidean plane, the concepts 


ta gcrange instability and complete instability are 
equivalent. 


PROOF Let this system be given by the differential 


equations 


OR 


5d e 
I 


g(x,y), 


ee 
1 


Enere f(x,y) and g(x,y) are continuous functions. 

Glearly if this system is completely unstable, then 
AEG is Lagrange unstable. So assume the system is Lagrange 
unstable and that there is a point p which is non-wandering. 
Then p € oe The peint p cannot bear critico ALO 
then cl (y’(p)) would be compact contrary to hypothesis. 
Likewise, there exists a neighborhood V of non-critical 
points about the point p. 


p 
Let p = a oP Since. p is not aili S a 
2 


either £(p, >»P,) #0 or 8(P, »Po) #0. Assume £(p, »Po) #0. 
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Since f(x,y) is continuous, there is a neighborhood W about 
p where f(x,y) X% O for all points a e W. Now consider 


the function H:R° + R defined by 


-] 
H(x,y) = tan (E Y. Clearly H(x,y) 1s continuous in 
3 0 


MN Since H(x,y) = Bann [2] = tan (SL), we may consider 
X 


the function H to be a "slope" function since it provides 
Meewith the slope of any trajectory at a given point. 

Now by the continuity of H(x,y) there is a neighborhood 
U' such that for all points E ín U', 
|IH(x,y) = H(p} »P,)| < e, where e > 0 is chosen such that 
U'C W. The relationship between the sets V, W, and U' 


WS shown in the following diagram. 


Y 
x 

Now since |H(x,y) - H(p, »P5) | < €e for all G) In Un 
we see that the neighborhood U' contains only points whose 
trajectories have slopes that differ from the slope of 
Y(p) at p by less than e. Now let & > 0 be so chosen 
meat 6 < € and 6 < 1/10. Now let 
U o) e U': |H(x,y) - H(p,.P5)| < 6}. The number 1/10 
was chosen so that we are assured that we can construct 


the following "transversal" through all the trajectories 


iene within U. 
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Let L be a line segment through the point p making 
an angle of 90% with the motion through p. Such a line 
can be constructed since we are working in the plane. Now 


every trajectory in U must cross this line L by our above 


menstruction of U. 


E 
A 


me claim that each trajectory crossing L in U must do 

so in the same direction as y(p). To see this, let re U 
and suppose y(r) crosses L in the opposite direction of 
y(p). Let p' be the point where y(r) meets the line L. 


Now let U-L = U U U igeese clear that Ul and U 


3° 1 5 are 
BS Joint sets. 





Since y(r) ana y(p) cross L in different directions, 
we may assume that for some small time t' > 0, p't' e€ U) 
amd pt' e U,. Now let L' be the segment of L which joins 
the points p and p' and consider m(L',t'). Since L' is a 


connected set, m(L',t') is a connected set. Also part 
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MDL"! will be in Us and part of L't! will be in U, 


Bene p't! e U and pt! e Un. 


ib 2 
IN 
Bu 
AU 
L 
Now L't' is a connected set so there is at least one 
Beene q such that q € L'N L't'. Thus it follows that 
q = qt' since it would be impossible that any other 
brajectory passing through L' at time t = 0 would cross 


the point q in time t = t'. For this to happen, we would 


Mae the following situation. 





Clearly this is an impossibility in U and hence we 
have q = qt'. But this implies that q is a periodic point 
and hence cl (y’(a)) is compact and we have a contradiction. 
Hence all trajectories cross L in the same direction. 

Now p € I*(p) so there exist sequences {x} in X 
and (Ge in R such that x, + p and t, + +» with xt NDE 
If any of the x's are periodic, then cl GG) is compact 
and we have a contradiction. Hence no Xn 1s periodic and, 


without loss of generality, we may assume each Xn iTeszon 
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L. Now for n sufficiently large, we may assume that 
et. e U since N. > p. Then it follows that there 
Faısts a Te R such that Ze) licctont io Cali a 
Point 2. 

It is clear that z lies between the points x and p, 


for if not we have the situation shown below and it is 


clear that in this case, x > DS an impossibrasiere 





Now, let L, be the segment of L which joins the points 


1 


x, and z. Since z = a + T), we consider the set 


{x [o,t, + TJU L,} = A, This is a Jordan curve and hence 
the interior of A together with A is a compact subset of 
ES. Let B denote this compact set. Now p £e B so we 
consider y'(p). We claim that pt e B for all te R. If 
this is not the case, then Ytp) must cross the set A 


+ 
at some point since yt) is a connected set. But y (p) 





cannot cross A anywhere along x_L0,t,_ + T] because different 
Meryeccories cannot intersect. Then the only alternative 


is that y'(p) must cross L 


1" 





Fut this too is impossible since to do so would mean 
Prossing L in the opposite direction from y(x). Hence 

e e B for all t e€ cae Hence v (p) c B and sincezbrTes 
compact, we have cl AE B. Therefore cl (y* (p)) 

is compact. Thus we have our final contradiction and the 
assumption that p € O) is untenable. Thus the system 

is completely unstable. We 
(The above proof was communicated to me through a private 


conversation with Dr. William Chewning.) 


We now continue with several characterizations of 


dispersive systems. 


t 


THEOREM 4.15 A dynamical system (X.R) is dispersiwe IE 


and only if Tx) = Y for each X € X. 


PROOF Let the system be dispersive. Assume there is a 
+ + 
point x e X such that J (x) #’ @ . Then ete e 


there are sequences lx, and ee such that x > X, 
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t, o and lS > y. This shows that for any neighborhoods 
U, of x and Uy Om ut N U # Ø as the element xt 

is contained in this intersection. But this contradicts 

the assumption that the system is dispersive because 9, is 
then positively recursive with respect to U,. Hence 

J” (x) = YM for all x e X. Conversely, let I” (x) = Y for 

all x e X. We claim that there exist neighborhoods Uy of 

x and U, of y and a T > 0 such that U t NA Do gto gel 

t > T. For if not there exist sequences {x} EN, es) 


such that X, >X, y + yand E, + TO, so that 


n me In 
y € Milz). Pursechis JS a contradiction. molmi larly Uy 
is not positively recursive with respect to U. Hence the 


system is dispersive. I4 


THEOREM 4.16 For a given (X,R,w) the following are equivalent. 
1. The space (X,R,7) is dispersive. 
MeO any two points x,y im A there ane neweaben— 
hoods U, of x and Uy of y and a constant Tez 70 
such that U, N Ut = Ø for all t, e ee 
any two distinct points x;y dn X, y £ ES 


+ 
a £ J (y). 
PROOF Let (X,R,T) be dispersive. Then for any two points 
x,y in X, there exist neighborhoods U, and Uy such that U, 


y 
We positively recursive with respect to U,- Hence for some 


is not positively recursive with respect to Do and U, is 


y 


t > T which implies EAS IATA uA Ut = 8 for 


> 0, U NÑ U,t = Ø forall t > T and utN U o a 


aait, |t| > T. 
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Assume statement 2 and suppose xX € O Then there 
exist sequences fy} and (t,) such that y, > y and t, > +o 
and Et, + x. For n sufficiently large we may assume that 
Di € Uy and fe € D and Ta O NS 
U, A Doi ahi ech is a contradiction. The proof that 
yz IT (x) is similar. 

Finally assume that for any two distinct points x,y in 
Y E lx) and x £ I’(y). Suppose for all neighborhoods 
BE and Dr and any T > 0 we have A pi A Y for some 
t > T. Since this is true for all neighborhoods, it will 
be true for a nested sequence of neighborhoods about the 
points x and y. Hence consider for each n, the neighborhoods 
S(x,1/n) and S(y,1/n). Now set T = T.. Then 
EDEA S(y,1)t, 7 § for some t, > T, Tees: 
Then there exists a to > T, + 2 such that 
B(x,1/2) f\ S(y,1/2)t, Z Ø. In this manner we have constructed 
sequences{x}, 2 Ka with x > x, ty? A T 
and Ya + y. We note that the process in eonstrüctingz.ne 
sequence tt} cannot terminate for if it did, then we have 
found a T, > 0 such that SELLO Nasen iron 
t > Ty contrary to our hypothesis. Thus we have A > y 
SO ye Fa) but this is a contradiction and hence Uy 
cannot be positively recursive with respect to Uy + A 


similar argument holds for Uy not positive recue 


with respect to U,. Hence (X,R,1) is dispersive. Ve 
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We now introduce the final concept in this section. 
It is that of parallelizable dynamical systems. We are 
headed toward one particular theorem and consequently many 
lemmas and theorems will be presented without proof since 
their proofs are Rosette ive in nature or are not particu- 
larly informative. However, the proofs may be found in 
Bhatia and Szego pages 48 through 53. Formally we begin 


with the definition of a parallelizable dynamical system. 


DEFINITION 4.17 A dynamical system (X,R,7™) is called 
parallelizable if there exists a set SC X and a homeo- 
morphism h: X > SxXR such that SR = X and h(xt) = (x,t) 


for every x £ S and t € R. 


To proceed in the study of parallelizable dynamical 


systems we need to develop a theory of sections. 


DEFINITION 4.18 A set SC X is called a section of (X,R,T) 
if for each x £ X there is a unique time t(x) such that 


Ox) € S. 


Merset S in Definition 4.17 is in fact a section of 
(X,R,7) for if z is any element in X, then z = xt for some 
xe Sandt eR. If this was not the case, then SR FX X. 
Then the value (-t) is precisely t(z) given in Definition 4.16. 
Īt is clear that zt(z) = z(-t) = x €e S. To see that (-t) 
is unique, assume there exists (-t') e R such that 


z(-t') = y e S. Then xt = yt'. Now using the homeomorphism 
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a a a PA 


h, we have h(xt) = h(yt') which implies Gt) A 
Thus we have x = y and t = t'. Thus t(z) is unique. 

Then to summarize, it ís clear that if a dynamical 
system is parallelizable it must have a subset S which 
contains a point of every trajectory so that SR = X. More- 
over, the set S can contain only one point of every trajec- 
Bory for otherwise +(x) would not be unique. To see that 
this last statement is true, suppose there is a trajectory 
wnich intersects the section twice. Then there are points 
Ma this trajectory such that x,y € S and xt = y for 
some t > 0. Then t(x) = 0 since xo = x € S but also 
me = y € S and hence t = t(x). Since t > 0, this contradicts 
the uniqueness of t(x). This implies in particular that 
moe dynamical system cannot have any critical points or 
periodic trajectories. To see how a parallelizable system 
compares with the other dispersive concepts defined in this 


section, consider the following example. 


EXAMPLE 4.19 Consider a dynamical system defined in Ro by 


the differential equations 


ax 


at = PY 
dy - 
ar 
where f(x,y) is continuous, and moreover f(x,y) = 0 whenever 


the point a? Sto Che lorn O with n a positive integer. 
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For simplicity we assume that f(x,y) > 0 for all other 


Points. This system is sketched in the diagram below. 








Now we form the system in which we are interested by 


deleting the sets 


x 
IL, = JE AA Y E ee es Dee 


from the plane rae This system is sketched below. 
Y 





€ A  _ __ __=—_——— 
i A A 
Y o 
It is clear that this system is dispersive since the 
only way two neighborhoods po of x and Uy ue. yy evel! lex= 
positively recursive with respect to each other would be 
if the points x and y lie on the same trajectory. Then 
clearly we could find a T > 0 for which un Ut = Ø or 
un Ur = ø for all |t| > T. Hence this system iss gene 
sive but it is not parallelizable since no section S cen 


exist for which t(x) would be continuous. This fact can 


be seen when it is realized that the only section possible 
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is one that looks like y = 1/x + c for any constant c 
together with some section which intersects the trajectories 
below the x-axis. If this "lower" section intersects the 
trajectory on the x-axis, then t(x) could not be continuous. 
If the lower section looks like y = -1/x + c, then t(x) 
would be continuous but no section will intersect the tra- 
jectory on the x-axis. Hence the system is notsperatlle lrzabre, 
We can get an example of a parallelizable dynamical 
system from the previous example by considering the 


differential equations 


ax 7s 
ae 
ay 
T 


This system is sketched below. 
N 


Clearly any line parallel to the y-axis is a section of 


the system. 


These examples also give an intuitive feeling for the 
notion of a parallelizable dynamical system in that all 


trajectories are in some sense parallel. 
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The function t(x) will be basic in our development. 
In general the function t(x) need not be continuous however 
for a dynamical system to be parallelizable there must 
exist a section where t(x) is continuous on X. This fact 


meeproved in the following theorem. 


THEOREM 4.20 A dynamical system (X,R,w) is parallelizable 


PROOF Suppose (X,R,T) has a section S with t(x) continuous 
on X. Then clearly SR = X. We now define the homeomorphism 
EXA * SXR by h(x) = (xt(x),-T(x)). To see that h is one- 
Zone let h(x) = h(y) for x and y in X. Then xt(x) = yt(y) 
and -t(x) = ~t(y). Hence xt(x) = yt(x) e S which implies 
Bat x = y for otherwise t(y) is not unique. Hence h is 
we oone. The function h is continuous by the continuity 
of t(x) and the continuity axiom. The inverse 
nd: sxR + X is given by nt (x,t) = xt and is clearly 
one-to-one and continuous. Thus h is a homeomorphism so. 
the definition of a parallelizable dynamical system is 
Satisfied. 

Now assume (X,R,T) is parallelizable. Then the set s 
in the definition is a section of X. Since for any X € X, 
x = yt for some y £ S and t e R, we set t(x) = -t. Then 


MET) = y c S. The continuity of I(x) follows 


from the continuity of h in the definition. // 
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To continue this development we need to define special 


pes of sections. This is done in the following definition. 


DEFINITION 4.21 An open set U in X will be called a tube 
if there exists a T > 0 and a subset SC U such that 

sI, E U and for each x £ U there is a unique t(x), 

It(x)| < tr, such that xt(x) e S. Here I, = Tr 1) 
also called a t-tube with section S, and S a (T - U)- 
Section of the tube U. If I, = R then U is an ®-tube and 


S an (e - U)-section. 


DEFINITION 4.22 Given an open *-tube U with a section S 

and t(x) continuous on U, and given sets N and K such that 
NCKCS where N is open in S and K is compact, we shall 
call KR the compactly based tube over K. Then indeed T(x) 


restricted to KR is continuous on KR. 


Proceeding toward our final result, the next theorems 


and lemma are presented without proof. 


THEOREM If X is locally compact and separable, and if 


nme hun 


every X € X is a wandering point, then there exists a 


ne A A ne 





d tubes 
countable covering {K,R} of X by compactly based tube 


KR each with t Ox) continuous on K R- 


LEMMA A compactly based o-tube U with section K of a 


dispersive dynamical system (X,R,") is elesea a A. 


We also need the following important topological 


theorem. 
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TIETZE EXTENSION THEOREM Let X be a normal topological 





SS Ą[Á ome cm Se ŻA ți ël 
e. 


space, let A be a closed subset of X and let f be a 


continuous function on A to the closed interval [a,b]. 





Then f has a continuous extension g which carries A into 


k.b]. 


The final lemma needed to prove our result is quite 


instructive and is completely presented below. 


LEMMA Let U; Us be two compactly based tubes of a disper- 


gee Se A A AAA -e ë D ug 
e m ë A AA — 


sive dynamical system with sections Ky > Ky and continuous 
functions t,(x) and 1,(x) respectively. If U, A Us * Y 


then U = U} U U, is a compactly based tube with a section 


A UA E E 


K DK and a continuous function t(x). Moreover, if the 


time distance between K, and K, along trajectories jam 


v A U, is less than 1 (> 0), the time distance between K 


and K, along trajectories in U is also less than T. 


PROOF In this proof, the following diagram willi be usera 


s UAU 
a 
U, fl 4 LLZ / 
J PZITISN ÍA — 
IR Ky AA 


nn 


By the previous lemma, Uy and U, are closed. They are 
invariant because each tube contains a section Ky and Ko 
respectively and K,RC U} and K,R E Us. Hence U N U, is 


invariant and closed. Further, Ko, A U, is compact and 


or 





non-empty, as shown in the diagram. Set S, = Ko A U, 
and S, = K, A U5. Any trajectory in uA U, intersects 
Sy wA actly one point and intersects 3, in exacte ly one 
Merny. Thus for any X e€ U, A Us, tT, (x) = 150%) + Ta (xT (X)). 
This is true because xT, (x) = XT5(X) (1, (XT 5 (x)= x(t, (X) + 
1], (xt5(x))) and there are no rest points or periodic 
trajectories in a dispersive system. The function Ty is 
Fonrtinuous on So which is compact and now we apply the 
Tietze extension theorem. This is possible since every 
metric space is normal. Thus the continuous Tunersıon Ta 
is extended to a continuous function t defined on Ko where 
t(x) = t, (x) for xe S,. If tT, (x) O fOr enre 
we have t(x) € (-1,T) for x € K,. 

Notice now that {xt(x): x € So} = S} and (x) being 
Memtinuous, {xt(x): xe Ko) is compact as K, is compact. 


Set K = KU (xix): xE Ko! and define t*(x) on 


KR = KR U KoR as follows: i 


t, (x) for Ne K,R 
T*(x) = 


To (X) + T(xto(x)) for x € K5R 


ı*(x) is continuous on KR and we need only verify chat 11 
m € un U,, then T} (x) = T,(x) + T(XT) (x)) which has 
already been proved since t(x) = T} (x) for x € S,. Thus 


the lemma is proved. A 
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Finally, we prove the theorem that is the main focus 


of this section on parallelizable dynamical systems. 


THEOREM 4.23 A dynamical system (X,R,7) on a locally 
momoact separable metric space X is parallelizable if and 
Emn if it is dispersive. 


PROOF Let the dynamical system (X,R,7) be dispersive. By 
Theorem 4.20 it is sufficient to prove that X has a section 
wen t(x) continuous on X. By a previous theorem, there 
is a countable covering UA) of X by compactly based tubes 
U with sections K, and continuous functions TDR We 
replace this covering by a like covering {u"} of compactly 
based tubes which we construct as follows. Set Ky =K, 
and U = u Beginning with u and Us we use the previous 
lemma to enlarge Kl to a compact set ce thus obtaining 

the compactly based tube ue = utu U, with (o) continuous 
on US. This leaves KY unaltered. Having found U" in the 
same manner, we take it together with UFI anda "conse lco 


> 
nl) continuous on y” L Now 


Di mitn K DK”, andr 
set S = OK, then X = SR and the function t(x) defined by 
T(x) = T (x) mor X E U? is continuous on X, with the property 
that xt(x) € S. Moreover t(x) is unique for each x € X. 

Thus X has a section S with continuous t(x) defined on X. 

The system (X,R,t) is thus parallelizable. 


Conversely, suppose the system (X Rom) is pemallelligable. 


Then there exists a section S with TÍ) continucús on A: 
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Suppose I) ZERO Some x EX. Tet y € Tx). 
Then since y e X, there exists y e S and ty e R such that 
yt = y. Now since J*(x) is invarianie O € TF) 
which implies that y e TF. ‘Then there exist sequences 


mix} and 20) with Zur, 


Y to, and Le = y. Now 


for each x_, there exists 2 e S and ta e R such that 
x = x. Similarly there exists xe S and u e R such 
Chat xt, = x. Then since x to > y, Jtesollews= cleat 
Ee t) > y, 
n ader for each fixed k, k= 1, 2, ..., the sequence 


{xt (t,)}. Since x,t, + xt,, we have by the continuity 


Axiom that xt (t,) > xt (t). Hence we may assume by 

eking subsequences if necessary that 

d(x_t (t,),xt,(t,)) < 1/k for k = 1, 2, ... Then it follows 
oo z 

that xt (t) + y since 


a(y,xt (t_)) < d(y,x,t (t,)) + a(x, t (t_),xt,(t 3) < 


acy ,x tt)» + 1/n. Now since xt, (t,) > y, let 


n 
A Dy + th Since E, Fee, it Lollowsechalzr 2252: 


Hence we have that_y e NR). 
A 





S ae nz 
Now for each n set le XT: Then yan sD SINCE Ta AS 


we may assume that To > n for each n. Let Uy be an 
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arbitrarily small neighborhood of y. Then since Ye TE 


3 


Y, E de for all n > N where N is a large enough integer. 





We know that T(x) is continuous and since (y) = 0, it 
follows that ya is Close to 0 for ali n= Ne Ene eRe 
may assume that [r(y < 1/n for each n > N. But this is 
impossible since El is unique and we have already 
determined that I) en for allen. Thus we Mare 
een radiction so Ít is impossible that y e UE and 
hence Mix) = @ for all x € X. Then by Theorem 4.15 we 


have that (X,R,T) is dispersive. Y 


The relationships between the various dispersive concepts 


is now conveniently summarized below. 





parolleliz able 






dispersive 






complete ly unstable 







Poisson unstable 





La qrange Unstab le 


a 





N= SPABT LIT DEREN 


In this section we will study stability theory first 
Brom the abstract point of view that we have established 
and then from an applied point of view that uses systems 
Soedifferential equations and studies the eigenvalues of 
such a system. To begin we need some basic definitions. 


moroughout this section, it will be assumed that the space 


Eee i ec eS ee 


X is locally compact and the set M is a non-empty compact 
Bubset of X unless explicitly stated otherwise. 


BERTNITTITON 5.1 With a given ME X we have 

l. the set AM) = {xe X: A(x) MA pie 

2. the set A(M) = {xe X: A (x) # Ø and A (x)c M}. 
The set AM) is called the region of weak attraction and 
the set A(M) is called the region of attraction of the set 
M. Moreover, any point in AM) or in A(M) is said tope 


weakly attracted or attracted to M respectively. 


From this definition it is clear that a point x is 
weakly attracted to a set M if and only if there is a 
sequence {tt} oo So + and d(xt, ,M) > 0. Also a 
meint x is attracted to a set M if and only Td AU O 
as G > to, 

The above attractor sets are important in the development 
of stability theory and the following theorem provides 


needed information about them. The following lemma is useful. 
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LEMMA 5.2 Let X be any metric Space and let x € X. Then 


Bor every te R, 


l. the ¡MES 


AT (x)t = AT(xt), 


set 
+ 
e set (a) = 3° (x)t = lee). 


N 
N 


PROOF The first equality in statements 1 and 2 follows 
from the invariance of the sets AT (x) and Fa) respectively. 


To see the second equalities, consider the case 


At (x)t AY (xt). Let 2 ¢ Moor. Maen there io anme A* (x) 
with z = yt and a sequence ice ave > aV 
Einen by the continuity axiom, xt, (t) > yt. But 

xt, (4) = xt(t,) and since t} > +, we must have that 

BL € A’ (xt). Thus AT(x)t C A (xt). The argument is clearly 
reversible so we also have that ATOxt) C AE Hence 


statement 1 is proved. The proof of the second inequality 


in statement 2 is entirely analogous. // 


EHEOREM 5.3 For any given M, A, (M) -> A(M) and the sets 


AM) and A(M) are invariant. 


ROOF Let x € AD. Then A(x) ž Ø and A (x) C M. Then 
clearly ATxx) N Mx BG. Hence x € ADE To see the 
invariance of the sets A, (M) and A(M), consider the case 
for A.M). Let x € A, (M) and let t e R. Now 

CKN Mx YB and from Lemma 5.72, A* (x) = AT (xt). Hence 
Next) N MX Ø which implies that xt € Am). Since t eR 
was arbitrary, we have that A, (M) is invariant. The proof 


of the invariance of A(M) is entirely analogous. 77, 
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Using these sets, we are now ready to define the various 


types of stability to be investigated in this section. 


DEFINITION 5.4 A given set M is said to be 
l. a weak attractor if A (M) is a neighborhood of M, 
2. an attractor if A(M) is a neighborhood of M, 
3. stable is every neighborhood U of M contains a 
positively invariant neighborhood V Of “ig 
4. asymptotically stable if it is stable and is 
an attractor, 


5. unstable if it is not stable. 


Geometrically, Mis stable if given any neighborhood U 
of M, there is a neighborhood Vve- uU such that any era eco) 
that is within V at time t = 0 remains within \ organ 
t> 0. If, in addition, all trajectories approach M as 
t increases, then M is asymptotically stable. Finally 
M is unstable if every neighborhood of M contains trajec- 
tories which are arbitrarily close to M at time t = O but 
which eventually leave the neighborhood for all time. These 


three cases are illustrated in the diagram below. 
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Nefnave seen specific instances of these cases in 
Several of our examples. For instance in the pendulum 
problem, (Example 2.1), the point (0) is stable; in the 
damped spring problem (Example 2.2), the point (5) alse 
asymptotically stable and in Example 3.2, if the system is 
meee led to the interior of the unit circle, the point 
P) is unstable. In this latter example, if x is any point 
Ne unit Circle, then the point x is a weak attractor 
mernough it is not stable. In order to see an example of 
Pamaceractor which is not stable, consider the following 
example: 

Let the planar dynamical system be defined by the 


differential equations (in polar coordinates): 


r(l - r), 


2 
ul 


sn a (6/2). 


D 
| 


\ 


ioe trajectories of this system are sketched below. They 
Asist of two critical points at (0) and (Ne a trajectory 
on the unit circle with (0) as the positive and negative 
ir set of all points on the unit circle. All other points 
Pave the point (e as their positive limit set and hence 

the point (2) is an attractor. It is not stable because 
every neighborhood will contain part ot pne undt ciTe 

where 6 > 0. Hence no positively invariant subset orans 


neighborhood can exist. 


TS 








We are now ready to prove some theorems about the 


concepts just introduced. 


AM OREM 5.5 If M is a weak attractor (respectively, an 





attractor), then the set AM) (respectively, A(M)) is an 


open neighborhood of M. 


PROOF Let N denote either of the sets A, (M1) or A(M). Then 
N is an invariant neighborhood of M. Consequently, bdy N 
Meeinvariant by Theorem 2.10, is disjoint from M, and is 
closed. Then for each x € N, MX) N MX Y, whereas for 
each . DS N(x)c bay N since bdy N is invariant 

and closed. Since bay NAVM = Ø, we conclude that 


Nf bdy N = @. Thus N is open. IA 


Mefnote that if a set M is stable then it is “the inver-— 
section of positively invariant neighborhoods and M is 
therefore positively invariant. Immediately from this we 


see that if the singleton {x} is stable, then {x} must be 
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Bositively invariant so conseqauently, x must be a critical 
point. 

he" next two theorems show different characterizations 
of stable sets but first we establish the following lemma 
mad its corollary which do not depend on the local compact- 


ness of X. 


LEMMA 5.6 Let X be an arbitrary metric space. Let x e X 


ME O. Then J (x)C J (w). 


+ 
PROOF Given w e A (x) and any y € I* (x), there exist 
sequences IA Sy es! > Wi, and oe and {x}, 
a Xs th to, cn > y. We may assume, if necessary by 
Smoosing subsequences, that tA -- a >in TOPE ACI 
Bonsider for each fixed k, k = 1, 2, ..., the sequence 
| 1 . 1 1 = 
Px tt By the continuity axiom x,t) >xt), k = 1, 2, 
We may, therefore, without loss of generality, assume that 
for each fixed k, A(x) x Ty) < 1/k for n > k. This shows 


1 1 1 1 ! < 
that x T) > w, because d(w,x T1)<d(w,xT,) A ee 


n n 


' ; vo (eet 2 
d(w,xt!) + 1/n. Now notice that xt, x T (ts u); 


1 1 I _ as : 
and xt. et,” w, and ti Ta n Hence 
y € en). As y € 3* (x) was arbitrary, we have 


"ic Tem, and the lemma is proved. 7% 
+ t ar 
COROLLARY Given M and X€ AM), then J (Yc T OG bag): 


PROOF Indeed for any x € A (M), S a MAS. Uet w 

be any element in (x) AN M. Then by the previous lemma, 

(a) E ou (ar). Bor De Jt (M) and we know that 

yO E DOD. Hence J’ (x) C J’ (M) C D (M). // 
i 





| 
| 


We now give an important charactenization o Ts taba y 


ei a set M. 
THEOREM 5.7 A set M is stable if and only if D*(m) = M. 


PROOF Let D’(m) = M and suppose if possible that not 
Hole. Then there is an e > 0, a sequence tx}; and a 
sequence {tt}; with t Oks d(x 5M) > 0. tana a(x t M) Re 
We may assume without loss of generality that £ > 0 has 

been chosen so small that S[M,e] and hence H(M,e) is compact 
Mere H(M,e) = {x e X: d(x,M) = e}. This is possible since 
X is locally compact. Further, we may assume that 


x, * X € M. We can now choose a sequence {t,}, O< tT <t 


n n 


Such that ne H(M e) n = I, 2; eee Since n Da 
compact, we may assume that Xani y e H(M,e). Then clearly 
ye D (x)C D (M), but y £ M. This contradiction shows that 
M is stable. Conversely, assume that M is stable. Then 
given any neighborhood U of M there is a positively invariant 
meiehborhood V of M wita Væ U. Sances tor anya e 

D(x) C cl WR’ for any neighborhood W of x, we get 

mx) Cc civ since V is positively invariant. § Thus 

DMC co: any neighborhood U ofr M MHence 

D’(M)C ech U: U ig a neighborhood of M} = M as M is 


+ 
compact. Since M ED N) always holds, we have D (M) = M. // 


The following remark shows what results when we combine 


several of the concepts under discussion. 
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an attractor and consequently asymptotically stable. 


PROOF Since M is a weak attractor, AM) is a neighborhood 
of M. We need to show that A, (M) E A(M). Let x e A 0). 
Then A (x) % 8 since A*(x)N M # @ and A*(x) © I*(x) c D (M) 
by the previous corollary. Since M is stable we have 
D’ (M) = M by Theorem 5.7 and consequently AT (x) € M. Thus 
x € A(M) and A(M) is a neighborhood of M. Thus M is an 
attractor. Since M is stable and is an attractor, M is 
asymptotically stable. 

The various implications that have been proved thus far 


are summarized in the following diagram. 


asympö ically 
stable 
Wear 
ater tiradas | 


We would now like to confine our attention to the study 










of stability for systems of differential equations. We 
shall begin by considering autonomous systems of constant 


coefficient linear equations; that is, systems of the form 


AX 


>d o 
n 


daly 





Mere A is an nxn real valued matrix. The origin in R” 
is a critical point for such a system, and when zero is 
not an eigenvalue for A (so that is when A is non-singular) 
En origin 1S the only critical point. The solutions of 
mes system can be described in terms of the eigenvalues 
Sten in the following manner. 

Let X = a + bi be a complex eigenvalue for the nxn 

real matrix A, and let E, be an eigenvector in 0% belonging 


A 
to A. Then the functions 


X (t) eĉ" (G, cos bt + H.sin bo, 


A 


X, (t) eĉ" (H, cos bt -fG sin Dnr 


À 


are linearly independent solutions of X = AX, where 





ERRE 
. Bee A 
G. = 5 and 
H, = 5 . (Kreider, Kuller, Ostberg, 
Damen 


Now when a + bi, b > 0, are complex conjugate eigenvalues 
of multiplicity m, this system has 2m linearly independent 


solutions constructed from functions of the form 


at,K 
e 


Gr t cos bt 


and 
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k 
H e?"t“sin bt, 


Mere k is an Integer such that O <k<n- l. Then 


assuming t > 0, 


at, k 


e 


2 


| [G e 


and since||G,|| is a positive constant, 
lim ||G eat Kos bt || SEO regen 
F700 A 


On the other hand | 1G, e**t*cos bt|| is unbounded if a > 0. 
eimilarly, 
lim | IH, e?"t*sin bie) Seo rear 
t> 
at 


but ||H,e®"t“sin bt|| is unbounded if a > 0. 


per cof complex conjugate eigenvalues for A will approach 
the origin of R” as t > +o 1f the real part of these eigen- 
palues is negative; it will depart arbitrarily far from 
Memorifin, if the real part of these eigenvalues is positive. 
Miche first case, the trajectories exhibit the properties 
mequired for asymptotic stability whereas in the second 
fase they exhibit instability. This reasoning clearly 


applies equally well to trajectories arising from real 
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elgenvalues for A. It remains to consider the case where 
A admits a pair of pure imaginary eigenvalues + bi. If 
these eigenvalues have multiplicity one, the corresponding 


solutions of the system are constructed from the functions 


G, cos DE and Hsin bt 
and the fact that cos bt and sin bt are bounded as 
t + to yet do not tend to zero implies that the resulting 
trajectories have the properties required for stability, 
but not asymptotic stability. On the other hand tne 
multiplicity of +bi is greater than one, the system has 
solutions involving functions of the form 

k 


G,t*cos bt and Ht sin poet 


with k > 1. When this happens, the origin is unstable since 
|t*cos bt| and It sin bt | 


are unbounded. Thus we have proved the following theorem. 


THEOREM 5.8 If X = AX is an nxn linear autonomous system 
bese coefficient matrix is nonsingular, then the origin 
WR is 

i asympcotically stable if the Tea o al 


or the eigenvalues of A are neri tiii; 
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2. stable, but not asymptotically stable, if A 


has at least one pair of pure imaginary eigen- 
values of multiplicity one, no "pure IMacinany 
eigenvalues of multiplicity excecdikik OMAN 228 
nomereenvalvues with positive real Dans 


3. unstable otherwise. 
Consider the following plane autonomous systems: 


im xX = y 2 >; 3. X= -X 


The coefficient matrix for the first system is 


© 
| 


p 
© 


with eigenvalues A = tl. Hence the origin is unstable 
since one of the eigenvalues has a positive real part. 


For the second system, 


This system is precisely that of Example 2.3 and has 
eigenvalues A = +i. Thus the origin is stable but not 


asymptotically stable. This is exactly as we would expect 
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My recalling the geometry of Example 2.3. 


Hor the third system, 


-1 0 
A = 
0 -1 
maas system has A = -l as its only eigenvalue and thus 


the origin is asymptotically stable. 


The above results can be used to obtain a complete 
description of the trajectories of a linear plane autonomous 


system of the form 


Di ® 
N 


ax + by 
y = cx + dy 
where the matrix 
a b 
= 
e q 


is nonsingular. This description depends upon the nature 
of the eigenvalues A > As of A which are the 7009. oz 


characteristic equation 
a Ea vera 
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Since the discriminant of this equation is given by 
A = (a + a) i 
a td)” = ‘(ad - beje= (a - 1d)" + Be, 


me follows that Az and A will be 

irreal and distinct if ie 0, 

real and equal It A = 0% 

complex eamjugates if A < 0. 
Moreover, when Az and A, are the roots of this system, We 
have (1 - AA - A>) = 0. Thus we have 


2 


 - CA, + Ay IA + ALA 


11 
O 
w 


I 2 


So 


ad =- bc = det A 


> 
N 


aus when A > 0, A, and Ao will have the same sign if and 


ji 
only if ad - be > 0, and will then be positive or negative 
mreorgemeras a + d is positive or negative. On the other 
Dand, if A < 0, Az = a+ Bi, 1) = a - BL BO 
atrd= 20. 

In this case A and Ay will be pure imaginary if and only 
if a+ d= 0. Otherwise they will have a nonzero real part 


which agrees in sign with a + d. 


Now, combining these observations with Theorem 5.8, we 


have the following description of the origin for a plane 
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autonomous linear system with constant coefficients. 


a tO 0 asymptotically stable 
MA O0: sa td> unstable 
atd=0 stable 


ad - be > 0 and a+d<0  asymptotically 
stable 
2. A > 0:€ ad = be > 0 and atd> o0 unstable 


ad - be < 0 unstable 


U) 
> 
il 


stable 
ad - be < 0 unstable 


i B = bc > 0 and a+dx< 0 asymptotically 
Thus far we have not discussed the cases where A is 
Singular. These cases are rather uninteresting but for 
completeness, they are presented here. 
Suppose A is the 2x2 zero matrix. Then x = 0 for all 
X € ES which implies that each point in Ro fa a critileal 
point and consequently each point is stable since every 
point within some neighborhood of a point, remains within 
nat neighborhood as t > +». 
Suppose A is a 2x2 matrix with zero as its only eigen- 


value. Then A must be of the form 


where a is some real number. 
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To find the eigenvectors Xo associated with A, we set 
AXo = 0. Then all eigenvectors are critical points and 


are given by the equation 


a -a X 
0 = lj = ax, - ax 
al 2 
a -a = 
2 
which implies that Xy = X. Thus all points on the Tine 


rough the origin with slope = 1, are critical points. 


Now 
Xx, = ax, = ax, 
X, = aX, -~ aX, 
so we have X = X, which implies that x, = Xy + c where c 


is any constant. Thus the remaining trajectories also have 
slope = 1, and intersect the X-axis depending onc. 
Their directions are shown in the sketch below. It is clear 


that all critical points are unstable. 
| ke 


ZZ 


en 





Thus far we have studied only linear systems. To 
study non-linear systems we need the forrowing idea or a 


linear approximation. Consider the system 
X = F(X) 


where X is a vector in R”. We shall show that whenever F 
mora function of class ct (that is, F has continuous first 
partial derivatives) in a region QR of R” containing the 


origin, then we may replace F by a function of the form 
JX + G(X), 


where J is an nxn matrix and G(X) is "small" in comparison 


with JX when ||x|| is small. Precisely, we state 


DEFINITION 5.9 Let F be as above and let X, be a point in 


O 
Meen a function L = L(X) is said to be a linear 


approximation to F at Xy if L is linear on R and 


[POD = FOR = LOD -l 


lim 
| [x-x, 11>0 IX - xl! 


Moreover, the existence and uniqueness of the linear 
1 
approximation can be guaranteed if F is of class C. 
[Ostberg, Kreider, Kuller pp. 417, 418]. This linear 


approximation is in fact the nxn Jacobian matrix of F at 
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every point Ao: specifically. 


oF ð F 


1 a 
IX, af IX, 
Ta: 
oF OF 
OX IX 
1 = 
X Xo 
where 
F (X) 
ECX) = 
F (X) 


and eR ss Oor all i. 


To gain a better understanding of the linear approxima- 


tion, we give some examples. 


We find the linear approximation about er for each of 


ere following systems: 


eer eX 


P. (x,y) 


P (x,y) 


oe 
N 


where Pz and P, are polynomials with P, (0,0) = P,(0,0) = 0 


Then P, (x,y) = P, (x,y) and F, (x,y) = P,(x,y). Now if 
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P. (x,y) = ax +t byte 


a 


2 
ee a AX + b,y° + higher-order terme, 


2 
>X E ey + higher-order terme: 


P,(x,y) = d}x + e,Y + fxy +4 


y 


we have 


a) 
Hx] 


= = a, + cy E ca9Xk + higher-order terms and 
oF, 

ao ES + ] = 

Sy dy c)x + eb y + higher-order terms 


Now evaluating these partials at the point Es we have 


JF JF 
Io ios 
ee eo igs En 


Similar results hold for the partials of Po. Thus we 


replace F(X) by JX + G(X) where 


2 2 
cy + + -+ 
ay Dy C} XY aX boy 
, and G(X) = 


2 2 
dá e £, Xy E dx + es) + 


Cy 
N 


Hence the linear approximation is given by the following 


system: 


JX. 


>d e 
il 
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2. X= sinx+e - 1 


y = Xy 


Now, sín x and e? have the following Taylor series 


expansions: 


3 5 7 
Sex =) = = + = - 77 + - 


y ye 


3 2 
+ ST + higher-order terms 


> 


Hence we have x = x + y - 


= 


and applying the procedure used in the previous example, we 


have that the approximating system is given by 


X = JX 
where 
1 l 
lees 
0 0 


To analyze this system as to stability, we note that the 
matrix J is singular and hence our theorem does not apply. 


Thus we must analyze J using the same method we used in 
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the two previous cases of singular matrices. We see that 
J has eigenvalues A = 0,1 and hence does not fit into the 
category of the other two examples. However, it is an 
easy task to determine that all eigenvectors lie on the 
fine Of slope -1 through the origin and are all critical 
Boints. from the approximating system, it is easily seen 
that y = Q and hence all trajectories not on the line of 
critical points are parallel to the x-axis. Their direc- 
tions are easily determined from the equations. They are 


shown in the following diagram. 


—— yee ee 
O nner YA 
un Zn 
a 


mt is clear that each critical point is unstable. 

One might question now whether the complete theory we 
developed for linear systems also applies to the approxima- 
ting systems for non-linear systems. Such is the case for 
asymptotically stable and unstable systems. Teaser. 
the linear approximating system is unstable, then the non- 
linear system is unstable. However, no definite conclusion 
can be made about stable systems. It is recalled that for 
a system to be stable, it could have no eigenvalues with 


positive real parts. A non-linear system may have 
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eigenvalues with small positive real parts which do not 
appear in its approximating system. These facts are 


Maeda” in the next section and are presented here Just 


for information. 
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VI. LIAPUNOV STABILITY 


The basic feature of stable theory a la Liapunov is 
that one seeks to characterize stability and asymptotic 
stability of a given set in terms of a non-negative real- 
valued function defined on a neighborhood of the given set. 
This technique was devised by the Russian mathematician 
A. A. Liapunov, and is known as Liapunov's direct or second 
method. It is based on the well-known fact that a physical 
system loses potential energy in a neighborhood of a point 
of stable equilibrium. More precisely, a point of stable 
equilibrium for a physical system is a point at which the 
potential energy of the system has a local minimum. This 
fact is known in physics as Lagrange's theorem. We proceed 


with the following lemma. 


LEMMA 6.1 Let the phase space X be arbitrary and let KC X. 


Let 6 be any continuous real-valued function defined on K 


A A A A A ie m m GE 
A GE | nn o l 


such that $(xt) < $d(x) whenever x[O,t]J c K, t > 0. Then 





Br for some x, cl (Y(x))C K, we have (y) = $(z) for 


Ewery y,Z € A’ (x). 


PROOF Now there are sequences tt} and Ku in R such 

meat t > +>, t. > +, and xt. > y, XT, > Z. We may assume 
n n n n 

by taking a subsequence that TÍ E, for each n. Then 


clearly o(xt_) = (xt) since xt, = xt (1, - ts 


a = Œ K. Thus proceeding to 
E > 0 and xt [0,T, tJ p g 
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Me limit we have by the continuity of 6, o(y) > ¢(z). 


A similar argument shows that $(z) > o(y). Thus the lemma 


ms proved. 77 


Me best known result on asymptotic stability is the 


following. 


THEOREM 6.2 A compact set MCX is asymptotically stable 
function E defined on a neighborhood N of M such that 
tee thx e M, we have Elx) = 0 and it xiae i ee 
have E(x) > 0, 
or x ¢ M, t > 0 and xí0,t1C N, we Have 


EAS ESO 


PROOF Assume that a function E as required iso ven. 

Choose a > O0 such that S[M,a] C N and is compact. Let 

m = min{E(x):x e H(M,a)} where H(M,a) = {x e X:d(M,x) = a}. 
Such a minimum exists because E is continuous and H(M,a) 

is a closed and bounded set. By statement 1 and the 
continuity of E we have m > 0. Set 

K = {x e S[M,a]:E(x) < m}. Then K is compact and because 

of statement 2, K is positively invariant. This establishes 
that M is stable as K is a positively invariant neighborhood 
of M. To see that M is an attractor, choose any compact 
positively invariant neighborhood K of M with KÆ N. Then 
or any xX € K, Go € K and since K is compact, 


+ 
cl (yrix)) eC K. Hence cl (y (x)) is compact and we have 
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Ø F Tc K, and Lemma 6.1 shows that E is constant on 

Br). But this shows by statement 2 tenat MG E Me Thús 

Mis an attractor and, consequently, asymptotically stable. 
Conversely, let M be asymptotically stable and A(M) 


its region of attraction. For each x ¢ A(M) define 
$(x) = sup{d(xt,M): t Sols. 


Indeed $(x) is defined for each x e A(M) because if 
d(x,M) = a, then there is a T > 0 with XLT, to) CE Na). 


This is because x e A(M). Thus 
E Sp a (xt TAO O 


Since d(xt,M) is a continuous function of t, é(x) is defined. 
Now d(x) has the following properties: $(x) = 0 for x e M, 
mx) > 0 for x Z M, and $(xt) < 6(x) for t > 0. This is 
clear when we remember that M is stable hence positively 
invariant, and that A(M) is invariant. Thus if (x) is 

Er ned for any x € A(M) it is defined for all xt with 

t e R We further claim that this 4(x) is continuous in 
my, Stability of M implies continuity of $(x) on M. 

For x ¢ M, this proof requires a knowledge of uniform 
sebo rs and it will not be presented here. However, 
this function is continuous in A(M) [Bhatia, Szego p. 67]. 


mies function may not be strictly decreasing along 
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trajectories in A(M) which are not in M and so may not 
satisfy statement 2 of the theorem. However we can obtain 


meescdesired function by setting 
E(x) = J o(xt)exp(-t) dt. 
0 


Clearly E(x) is continuous and satisfies statement 1 in 
A(M). To see that E(x) satisfies statement 2, let x £M and 
t > 0. Then E(xt) < E(x) holds because p(xt) < o(x) holds. 
To rule out E(xt) = E(x), observe that in this case we 


must have ¢(x(t + t)) = ¢(xt) for all t > 0. Thus in 


I 


emticeular, letting T 0, t, 2t, ... we get de AO 
met, 2, ... But asymptotic stability of M implies that 

for x € A(M), d(xt,M) > O as t > +œ. Thus 

¿(x(nt)) > 0 as n > o, as 4(x) is continuous. This shows 
that d(x) = 0. But as x £ M, we must have d(x) > 0 and 


we have a contradiction. Hence E(xt) < E(x) for x £ M 


and SOS Tus statements 1 and 2 are satisfied. // 


This theorem en about the size of the region 
w traction of M. Thus if a function E(x) is known to 
rst in a neighborhood N of M, we need not have either 
NCA(M) or A(M) < N. This means that the above theorem 
cannot immediately be stated as a theorem on global 
asymptotic stability. The following example helps to 


Pmelwstrate this point. 
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EXAMPLE 6.3 Consider a dynamical system defined in the 


n plane by the differential equations 


cea EXSY De, 
DI Y 
where 
x is xy? ze 
f(x,y) = 


2x y ~ x if sy < 


The trajectories of this system are sketched below. 





“ye 





Ane origin (0) is asymptotically stable, with the set 


(C8) exey® <1} as its region of attraction. Consider now 


he function 


SS 








Gy ee 


This function satisfies the conditions of Theorem 6.2 in 


the whole plane. To see this we note 


A ros F _ 3E 
at at + ae o ae. en =Y F 


It can be verified that E) <0 for ali (2) 7 O), which 
implies that E satisfies the conditions of the theorem in 
every neighborhood of the origin. But not every neighbor- 
hood of the origin is contained in the region ob attraction, 
nor does it contain the region of attraction. 

In view of this example, we wish to establish conditions 
necessary to ensure global asymptotic stability. This is 


done by starting with the following definition. 


DEFINITION 6.4 A continuous real-valued function E(x) 
defined on a set N C X will be said to be uniformly unbounded 
on N if given any a > 0 there is a compact set KCN, 


Mazen, such that E(x) > a for x £ K. 


The following three theorems will be stated without 
proof since the proofs are very laborious. The proofs may 


around in Bhatia and Szego, pp. 70 = 73. 


THEOREM 6.5 Let MCX be a compact asymptotically stable 


mm ne 
= =e = 


Eemeanad let M be invariant. Then there exists a continuous 


uniformly unbounded function E(x) defined on A(M) such 


that 
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l. for x e M, E(x) = O and for x e T, Dw > 20% 


2. for all x e A(M) and t € R, E(xt) = e7 E(x). 


BEEOREM 6.6 IEMCKX is any compact asymptotically stable 
set, there exists a continuous, uniformly unbounded function 
E(x) on A(M) such that 

l. for x eM, E(x) = O and for x £ M, E(x) > O, 


eo mame x EM and t > O,nElxt)ee Bes 


continuous uniformly unbounded function E(x) defined on 
an open neighborhood N of M such that 
Uor x e M, Elx) = 0 and- for xM, ECO an, 
Or? x £ M, t > 0 and x0 tre TAE E 
Then M is asymptotically stable and NC A(M). If in 
addition, any condition guaranteeing the invariance of N 


holds, then N = A(M). 


These three theorems are very Similar in their conclu- 
sions but we can see that the hypotheses are becoming more 
general. The next results on global asymptotic stability 


follow from these theorems. 


THEOREM 6.8 A compact invariant set M C X is globally 
asymptotically stable if and only if there exists a continu- 
ous uniformly unbounded function E(x) defined on X such that 


l. for xeM, E(x) = 0 and for x £ M, E(x) > 0, 


=t 
2, for all x e X and t e R, E(xt) = e E(x). 
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PROOF The sufficiency follows from Theorem 6.7 since X is 
f@eeanvariant neighborhood of M. The necessity follows from 


Theorem 6.5. I7 


THEOREM 6.9 A compact set MC X is globally asymptotically 


Mibounded function E(x) defined on X such that 
or ze, EX) = O and for x Z M, E(x) > 0, 


a M and € > 0, E(xt) < E(x). 


PROOF The sufficiency follows from Theorem 6.7 and the 


necessity follows from Theorem 6.6. // 


These theorems allow us to determine whether a set is 
a vme otically stable or not without having to solve the 
System. However, if a set is asymptotically stable, it 
may not be a trivial problem to display the continuous 
uniformly unbounded Liapunov function. As in the last 
section, we now focus our attention directly on systems 
of differential equations in order to apply some of the 
above theory. 

Thus we return to the general case of an autonomous 


system 
X = F(X) 


defined in a region Q of R with a critical point at Xo 


in 2. Without loss of generality, we may assume that the 
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Pe 


priiitcal peint is he origin 0 of R? for if it i ANOT, a 


Simple change of variable X* X - Xo will translate the 


wein to the point Xo without changing the form of the 


system. 


Again we seek a real-valued function E(X) of class ct 


defined on some region 2 of the origin such that E(X) > 0 


Pe 


and E(X) = 0 if and only if X = 0. We will refer to E as 
the energy function because of our previous physical 


interpretation. If X = X(X_ st) iS sUMemsolut Lon or 


X = F(X) 


X(0) = Xo » 


the time rate of change of energy, cE along the trajectory 


ained by this solution is 


JE 
— F 
19%, al 


VE + X = VE-F = 


tims 


al 


Thus if this en is negative in 2, a particle moving 
along any trajectory defined by the system will be dissipa- 
ting energy as it intersects the surface E(X) = e. Conse- 

wW aay a particle which enters the region enclosed by 

such” a súrface can never gather enough energy to escape, 

and its trajectory must remain in that region. This implies 
taauethe Origin will be asymptotically stable, for the 


given system. 
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These remarks motivate the following definition. 


DEFINITION 6.10 Let Q be a region of R? containing the 
origin, and let E = E(X) be a real-valued function of 
class cs in 2. Then E is said to be positive definite if 
Iran R, E(X) > 0 and 
2, the equality E(X) = 0 holds if and only if 


ag 


X= 0. 


n 
If in addition we have VE-F= p 22 Fr < 0 everywhere 


cee 
mame taen E 1s said to be the Liapunov function for the 


autonomous system X = F(X). 


In terms of this definition, we can now form some 


conclusions about the stability of the origin. 


Hoy 


x = F(X), 


~~ 


F(X) a O 


provided there exists a Liapunov function E for the system. 


PROOF Assume the existence of the Liapunov function for 
the system and consider the set 

s[0,r] = {Xe R": a(0,X) <r}. Let ropemso chosenmular 
s[0,r] lies in Q. Since s[0,r] is compact and E is 


continuous and positive definite in Q, it follows that E 
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achieves its minimum value on s[0,r]. Let m be this 
minimum value. Clearivem iS positive. 

Again since E is continuous, there exists an e > O such 
that |]E(x) - E(0) | | < m whenever ||X - o| | < e, or simply 
stated, E(X) < m whenever ||X]| < e. 

Consider the set o Let Xo 
be any point in S(0,€) aná let 
X(X,,t) be the unique solution 
of x = F(X) whose trajectory 


passes through Xo at time 





t = 0, 

Now we claim that the trajectory 
must remain within S[0,r] for all t > 0. Assume that this 
is not the case. Then ||X(X,,t,)|| = r for some t,> 0, 
Since the trajectory is a connected set. Thus 
E(X(X),t,)) > m but this implies ig is increasing with 

JE 


time and that — > 0. Now —= — ——= VE-F. 


yt SE aX, dE Hence 


VE-F > 0. But this contradicts the assumption that 
NOAA, since this inequality implies that E is 

a non-increasing eunetionzoer t alóne theatra cCtory in 
question. Thus the trajectory remains within S[0,r] for 
all t > O and hence the origin is stable since s[0,r] is 


a positively invariant neighborhood of the origin. 7 


This argument suggests that if E actually decreases 


along trajectories in a neighborhood of the origin, then 
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the origin would be asymptotically stable. The next theorem 


mses this fact. 


THEOREM 6.12 If E is a Liapunov function for X = F(X) with 
miemeroperty that -VE-F is positive definite in @, then 


the origin is asymptotically stable. 


PROOF We have seen that there exists an e > 0 with the 
property that E is a non-increasing function along any 
trajectory X(X),t) with | [xoll <e. Since E is positive 
definite in £, we claim that E approaches a non-negative 
A Eg fome nus trajectory as t 7 œl To see this, we 
form the sequence {E(X(X_,t,,))3 where t, > @. 

Now E(X(X) ty)? assumes values in [0 ,E(X,)]. Since 
[0,E(X))] is a compact set, every sequence in [0,E(X,)] 
has a subsequence which converges to a point E, in [0,E(X,)]. 
Moreover E is non-negative. 

We will be done if we show that the limit is zero when 
-VE-F is positive definite in 2. Then the fact that E is 
zero only at the origin will imply that the trajectory 
approaches the origin with increasing time. 

Mute Contrary, suppose Eo > 0. Then since E is 
continuous there exists a ô > O such 
that |]E(x) - E(0)|| < E, when 
MEAN consider the set 
Ios]. Now the trajectory 
X(Xq >t) cannot enter the sphere 


S(6,6] for if it did, then 
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E(X(X¿>8,)) < E, for some k > 0. But since E is asnon- 
increasing function, it could never leave S[0,6] and thus 
could not approach Eo as a limit. 

since is continuous and E is of class c+ mE, 
Ns Continuous and hence assumes its minimum value on 
the "annular" region 8 < ||X||] < r where r is as in the 
meeor Of the previous theorem. This fact is true because 
this "annular" region is compact. Let m be this minimum 


value of -VE-F. Since -VE-F is positive definite, it 


Ml lows that m > 0. Thus 


dE(X(X_,t)) 


TE IA lA 


Now 


E JE(X(X, ,7)) 


T OT 


EUA td) EQU) = J 


by the fundamental theorem of the calculus. But 


E JE(X(X¿>T)) t 
A at < MURANO =" Mb. 
0 i 0 


Thus we have E(X(X,,t)) < E(X) zmt. 
Since the right hand side approaches -~ as t + », this 


implies that E(X(X,) ,t)) < Omfor ele, where 


EO) 
ty = 0% put this contradicts the fact that E is positive 





m 
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ference inf. Thus Eo = 0 and the origin is asymptotically 


stable. Vy 


Thus far we have presented theorems only on stability 
SieasympcOovic Stability of the origin. However it is also 
important to know under what conditions the origin is 
Maoevaowle., Consequently we present the following instability 


theorem known as Liapunov's First Instability Theorem. 


i 


in 2 and suppose that E(0) = 0, E is positive definite and 
Bees able to assume positive values arbitrarily near the 


origin. Then the origin is unstable. 


PROOF Let S[0,r] and s[0,€] be spheres about the origin 
such that e < r and S[d,r]c 2. Since E(X) and the first 
partials of E(X) are continuous in N, we have that E(X) 

is bounded in S[0,r]. Hence E(X) < k for all X e S[Ū,r]. 
Choose any Xo € s[0,e] such that E(X? > 0 and let X(X,,€) 
be the unique solution of x = F(X) whose trajectory passes 


through X, at time t = 0. 


O 
Now X(Xp 8) must eventually cross bdy s[0,r] ior af eae 
then E(X(X) ,t)) approaches a limit p in [LE(X,),k]. Since 
E(X(X,>8)) + k as t +», it follows that Be + 0. 
Now E(X) = VE-F and since E and F are continuous, we have 
that E(X) continuous. Then since E(X(X, +4) en ee 
follows that X(X 95%) + 0 since E is positive definite. 


Pais bs a contradiction Since if X(Xp€) + 0 ast +o, 
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must follow that E(X(X,>t)) + 0 because E is continuous. 
Thus since E(X,) > 0 we have E(X(X, >8)) OU DUG hs aS 
a contradiction. Hence X(Xp st) must eventually cross the 


boundary of S[0,r] and thus the origin is unstable. Vy 


The important feature of these theorems is that they 
enable us to determine the stability or instability of 
Seeulcal points without actually solving the system. This 
alt icularly important in cases where the solutions are 
moss ible to Obtain in closed form or are difficult to 
analyze. We can still gain valuable information on their 
mooi yy provided that we can construct a Liapunov function 
for the system. Although there is no general method for 
constructing these functions, we shall see that this is 
fairly easy to do in many cases. To solidify our under- 


standing of this notion, we consider the following example. 


EXAMPLE The function E(x,y) = x° + y* is a Liapunov 


function for the system 


ee Bert xy, 


-y + xy“. 


< 
N 


Clearly, E is positive definite and of class o in the 
entire x,y-plane. Moreover, since 


TOE JE _ 
VE = zx êl + ne Ce exe, + eye, and 
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F(x,y) = (-x + xy)e, si xy“ )e, 
we have 
VE-F = _2(x° + yá) + 2x y + oxy? = 2(x° + ae er 


Hence -VE-F is positive definite in the region xy < 1, and 
E satisfies all the requirements of Theorem 6.12 in this 
feerom. thus the origin is an asymptotically stable 


em ical point for this system. 


ne shall mow describe a method of constructing a 
Liapunov function. To do so we begin by constructing a 
Liapunov function for the constant coefficient linear 


autonomous system 


>d o 
ii 


AX 


when all onee icenya Mues Or A have negative weal parts. 
It is noted that we previously determined that the origin 
In ano obmecall stable for this system. Thus we will 
Doe aan anything new from the construction but it will 
io sito extend some of our earlier results to non-linear 
systems. The construction goes as follows. 

Let Cys tees = denote the standard basis vectors in 


> ao cach 1, 1 < i <n, let 


R 
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X(t) 


X, (t) E X, Ce, st) = 


(kt) 


Xin 


denote the solution of this system which satisfies the 


mortal condition X, (0) = ej: Then if 
Y = y€} a T YEs is any vector in E, the fume wom 
X(Y,t) = yX (t) et Ya) 


ehe solution of the system which satisfies the initial 
condition X(0) = Y. We now set 


oo 


HCY) = J IIX, t) |1 dt: 
0 


and note that if this integral converges, then.E(Y) is 
Posmemvedetinive. It is clear that E(Y) > 0. Also if 
Y = 0, then for all a <n. We haye ea O and hence 
x(0,t) = 0. Thus E(0) = 0. Now suppose E(0) = 0. This 
implies that X(Y,t) = 0. Hence y, = 0 for all i,1<icn 
since otherwise X, (t) must be zero for some i. But then 
X, (0) # ey: Thus we must have that Y = 0. 

inemproor chat E(Y) is defined for all Y in R? can be 
found in Ostberg, Kreider, Kuller, p. 415. 

Finally we show that E is in fact a Liapunov function 


for this system. We compute the value of E(X(X,,8)) where 
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X(X,,€) foe tnemsOlUbion OF X ="AX saisi ying X(X, 50) = Xp: 


Hence we have 


oO oo 


7 a e BCE t)||° as 


E(X(X,»€)) 


00 


S | 1x0) 115 du. 
E 


11 


Hence 


gE(X(X,>t)) ad 


JE at 4 ea a 


as required since VE:*F = 9E/3t along the trajectories of 
X= AX. Moreover, since the real-valued function xp] 
is positive definite on R”, E satisfies all the hypotheses 
of Theorem 6.12 and again we have proved that the origin 
is asymptotically stable. 


We now use this method to construct a Liapunov function 


for the following linear system: 


y = -2y. 


This system has the coefficient matrix 
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with eigenvalues -1, -2. In this case the solution X(t) 


and X, (t) such chat x, (0) = €); X, (0) = e, are given by 


o 0 


Ne ee x. (u = 
il , 2 
0 pet 


enus if Y = xe, + ye,; 


00 
E(Y) = / | [xe "e + ye“"e,||° at 
> ISA N dt 
0 
ne yo 
A 


The fact that this function satisfies the hypotheses of 
Theorem 6.12 is easily verified and thus the origin is 


eer movoultcally stable. 


We now wish to extend these results to non-linear 
A stems. Yor this, we recall the notion of the linear 
approximation which was developed in the last section. 


We recall that in the system 
X = F(X), 


the function F may be replaced by a function of the form 
JX + G(X), where J is the nxn Jacobian matrix of F and 


G(X) is small in comparison with JX when ||X|| is small. 
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Now let F be a function of class gt in a region Q of 
R” containing the origin, and suppose that P (0) = 0 and 
the real parts of all of the eigenvalues of the Jacobian 
Mapa J of F evaluated at the origin are negative. We 


NOw show that under these conditions the Liapunov function 


E = E(X) for the constant coefficient linear system 


X 


JX that has been previously constructed is also a 


Liapunov function for the non-linear system 


X = F(X). 
Now by our earlier results, we know that VE-(JX) < e 
Hence 

VE(X)-F(X) = VE(X)-LIX + G(X)] 


= VE(X)-JX + VE(X)-G(X) 


-11x11% + |]veGo!1 110001]. 


[A 


This last statement is a consequence of the Cauchy-Schwarz 
inequality. Now we also know that ||VE(X)|]| < 2x]|X]], 
since VE is a linear transformation on a finite-dimensional 


Space. Thus we have 


VE(X) F(X) < -[ [x1]? + 2k} 1x1] 1100011. 
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We now use the fact from the definition of a linear 


E(X) = JX 


approximation that lim y 


|Ix||>o 


follows that lim G(X) = 0. Hence we can find a 


E 


OR Bab cen it 


sphere S[0,6] about the origin in R? such that 


IL 
eco II < q 1x1 
Mor all X e Some Then within this sphere we have 
a ANA ll El alo = 311116 
= : Uk 2 


and indeed E is a Liapunov function for the system 

x = F(X). Moreover, since this inequality also implies 
that -VE-F is positive definite in S[0,6], the hypotheses 
of Theorem 6.12 are satisfied and we have established the 


following result. 


THEOREM 6.14 Let F = F(X) be of class CÌ in a region of 
R? containing the origin, and suppose that F(0) = 0. Then 
the origin is asymptotically stable for the system X = F(X) 


meee SS ee 


Tecwmalvabedsan the Origin have negative real parts. 


As before, there is a companion result when all of the 
elgenvalues of the Jacobian matrix of F have positive real 
Patos aia that case the origin is a point of complete 


instability for X = F(X). 
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ÉS statements verify the comments in Section V 
about the conclusions that may be drawn from a linear 
approximation. Finally we show by means of an example 
that no definite conclusions can be made from a system 
whose linear approximation shows that the origin is a 


Amable critical point. 


EXAMPLE Consider the plane autonomous system 


y =- xf(x,y) 


> oo 
It 


-x - yf(x,y) 


do. 
N 


where f has a convergent power-series expansion in a 
neighborhood of the origin in Ro andere 05220. 


Taking the linear approximation, we have 


JX 


Ps 
i 


where 


© 
| 


In this case J has eigenvalues +i so by Theorem 5.8, the 
wena s stable for the system X = JX. 


Now let E(x,y) = ae ne Then 
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wr = Go (PAER) E 20 Dran. 
Hemee for this system the origin is 
1. stable when f(x,y) > 0 in some neighborhood 
of (0); 
2. asymptotically stable when f(x,y) is positive 
definite in a neighborhood of (0); 
3. unstable when f(x,y) < 0 in every neighborhood 


OR oon 


tnis conpletes our Study of Liapunov functions and 
stability theory in general. At this point, the reader 
should have sufficient background in this theory to be 
neo read and Understand further results in this area. 
Specifically recommended for further reading is Chapter 
VITI, cl -Liapunov Functions for Ordinary Differential 


Equations, in Bhatia and Szego. 
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